Hyalomma anatolicum, an Anatolian hard tick is a well-recognized vector involved in the transmission of various pathogens to animals and humans. The present study elucidated the population structure and haplotype network of H. anatolicum based on the mitochondrial large subunit ribosomal RNA (16S rRNA) gene sequence. The population structure and haplotype network analysis of 75 sequences archived in the GenBank, including the 15 sequences generated herein, yielded 24 haplotypes. Haplotype 1 (Hap_1) was the predominant haplotype consisting of 45 sequences from India, China, Pakistan, Turkey, Egypt, Iraq, and Tajikistan. The complete haplotype network exhibited a stellate conformation, highlighting a recent population expansion. The overall dataset, together with the sequences corresponding to India, China, and Pakistan, showed a high haplotype (0.638 ± 0.065, 0.671 ± 0.103, 0.753 ± 0.099, and 0.854 ± 0.061, respectively) and low nucleotide (0.00407 ± 0.00090, 0.00525 ± 0.00196, 0.00680 ± 0.00233, and 0.00453 ± 0.00056, respectively) diversity, further emphasized a recent population expansion. The neutrality indices including Tajima\'s D, Fu and Li\'s D, and Fu and Li\'s F for the complete dataset (- 2.661, - 6.008, and - 5.649, respectively) as well as for the sequences from India (- 2.223, - 3.414, and - 3.567, respectively) were negative, suggesting deviation from neutrality and a recent population expansion. The present study provided novel insights into the population structure and haplotype networks of H. anatolicum based on the mitochondrial 16S rRNA gene, and the different tests inferred a low genetic differentiation and suggested a recent population expansion of this economically important tick species.

RNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community\'s ability to interrogate both forms of RNA variation simultaneously.
We employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis.
Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

BACKGROUND: The cashmere goat industry is one of the main pillars of animal husbandry in Inner Mongolia Autonomous Region, and plays an irreplaceable role in local economic development. With the change in feeding methods and environment, the cashmere produced by Inner Mongolia cashmere goats shows a tendency of coarser, and the cashmere yield can not meet the consumption demand of people. However, the genetic basis behind these changes is not fully understood. We measured cashmere traits, including cashmere yield (CY), cashmere diameter (CD), cashmere thickness (CT), and fleece length (FL) traits for four consecutive years, and utilized Genome-wide association study of four cashmere traits in Inner Mongolia cashmere goats was carried out using new genomics tools to infer genomic regions and functional loci associated with cashmere traits and to construct haplotypes that significantly affect cashmere traits.
RESULTS: We estimated the genetic parameters of cashmere traits in Inner Mongolia cashmere goats. The heritability of cashmere yield, cashmere diameter, and fleece length traits of Inner Mongolia cashmere goats were 0.229, 0.359, and 0.250, which belonged to the medium heritability traits (0.2 ~ 0.4). The cashmere thickness trait has a low heritability of 0.053. We detected 151 genome-wide significantly associated SNPs with four cashmere traits on different chromosomes, which were very close to the chromosomes of 392 genes (located within the gene or within ± 500 kb). Notch3, BMPR1B, and CCNA2 have direct functional associations with fibroblasts and follicle stem cells, which play important roles in hair follicle growth and development. Based on GO functional annotation and KEGG enrichment analysis, potential candidate genes were associated with pathways of hair follicle genesis and development (Notch, P13K-Akt, TGF-beta, Cell cycle, Wnt, MAPK). We calculated the effective allele number of the Inner Mongolia cashmere goat population to be 1.109-1.998, the dominant genotypes of most SNPs were wild-type, the polymorphic information content of 57 SNPs were low polymorphism (0 < PIC < 0.25), and the polymorphic information content of 79 SNPs were moderate polymorphism (0.25 < PIC < 0.50). We analyzed the association of SNPs with phenotypes and found that the homozygous mutant type of SNP1 and SNP3 was associated with the highest cashmere yield, the heterozygous mutant type of SNP30 was associated with the lowest cashmere thickness, the wild type of SNP76, SNP77, SNP78, SNP80, and SNP81 was associated with the highest cashmere thickness, and the wild type type of SNP137 was associated with the highest fleece length. 21 haplotype blocks and 68 haplotype combinations were constructed. Haplotypes A2A2, B2B2, C2C2, and D4D4 were associated with increased cashmere yield, haplotypes E2E2, F1F1, G5G5, and G1G5 were associated with decreased cashmere fineness, haplotypes H2H2 was associated with increased cashmere thickness, haplotypes I1I1, I1I2, J1J4, L5L3, N3N2, N3N3, O2O1, P2P2, and Q3Q3 were associated with increased cashmere length. We verified the polymorphism of 8 SNPs by KASP, and found that chr7_g.102631194A > G, chr10_g.82715068 T > C, chr1_g.124483769C > T, chr24_g.12811352C > T, chr6_g.114111249A > G, and chr6_g.115606026 T > C were significantly genotyped in verified populations (P < 0.05).
CONCLUSIONS: In conclusion, the genetic effect of single SNP on phenotypes is small, and SNPs are more inclined to be inherited as a whole. By constructing haplotypes from SNPs that are significantly associated with cashmere traits, it will help to reveal the complex and potential causal variations in cashmere traits of Inner Mongolia cashmere goats. This will be a valuable resource for genomics and breeding of the cashmere goat.

Since the COVID-19 pandemic, the diversity of clinical manifestations in patients has been a tremendous challenge. It seems that genetic variations, as one of the players, contribute to the variety of symptoms. Genome-wide association studies have demonstrated the influence of certain genomic regions on the disease prognosis. Particularly, a haplotype at 3p21.31 locus, inherited from Neanderthals, showed an association with COVID-19 severity. Despite several studies regarding this haplotype, some key variants are not sufficiently addressed. In the present study, we investigated the association of rs17713054 at 3p21.31 with COVID-19 severity. We analyzed the genotype of 251 Iranian COVID-19 patients (151 patients with asymptomatic to mild form as control and 100 patients with severe to critical symptoms without any comorbidities as case group) using the ARMS-PCR method. Results demonstrated that the A allele confers an almost twofold increased risk for COVID-19 severity (P value = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (P value = 0.013). In conclusion, the A allele in rs17713054 was a risk allele in Iranian patients and was independently associated with COVID-19 severity. More studies are beneficial to confirm these findings in other populations and to develop strategies for risk assessment, prevention, and personalized medicine.

OBJECTIVE: To investigate the origin of Biomphalaria straminea in China, so as to provide insights into assessment of schistosomiasis mansoni transmission risk and B. straminea control.
METHODS: Guanlan River, Dasha River, Shenzhen Reservoir, upper and lower reaches of Kuiyong River, and Xinzhen River in Shenzhen, China, were selected as sampling sites. Ten Biomphalaria samples were collected from each site, and genomic DNA was extracted from Biomphalaria samples. DNA samples were obtained from 15 B. straminea sampled from 5 sampling sites in Minas Gerais State, Pará State, Federal District, Pernambuco State, and Sao Paulo State in Brazil, South America. Cytochrome c oxidase I (COI) and mitochondrial 16S ribosomal RNA (16S rRNA) genes were sampled using the above DNA templates, and the amplified products were sequenced. The COI and 16S rRNA gene sequences were downloaded from GenBank, and the sampling sites were acquired. All COI and 16S rRNA gene sequences were aligned and evolutionary trees of B. straminea were created based on COI and 16S rRNA gene sequences to identify the genetic similarity and evolutionary relationship between B. straminea samples from China and South America.
RESULTS: A total of 60 COI gene sequences with a length of 529 bp and 3 haplotypes were obtained from B. straminea sampled from China. There were 165 COI gene sequences of B. straminea retrieved from GenBank, and following alignment with the above 60 gene sequences, a total of 33 haplotypes were obtained. Phylogenetic analysis showed that the three haplotypes of B. straminea from China were clustered into one clade, among which the haplotype China11 and three B. straminea samples from Brazil retrieved from GenBank belonged to the same haplotype. Geographical evolution analysis showed that the B. straminea samples from three sampling sites along eastern coasts of Brazil had the same haplotype with China11, and B. straminea samples from other two sampling sites were closely, genetically related to China11. A total of 60 16S rDNA gene sequences with approximately 322 bp in length were amplified from B. straminea in China, with 2 haplotypes identified. A total of 70 16S rDNA gene sequences of B. straminea were captured from GenBank. Phylogenetic analysis showed that Biomphalaria snails collected from China were clustered into a clade, and the haplotype China64 and the haplotype 229BS from Brazil shared the same haplotype. The 49 16S rDNA gene sequences of B. straminea from 25 sampling sites in southern Brazil, which were captured from GenBank, were included in the present analysis, and the B. straminea from 3 sampling sites shared the same haplotype with China64 in China. Geographical evolution analysis based on COI and 16S rRNA gene sequences showed that B. straminea sampled from eastern coastal areas of Brazil shared the same haplotypes in two gene fragment sequences with Biomphalaria snails collected from China.
CONCLUSIONS: The Biomphalaria snails in China are characterized as B. straminea, which have a low genetic diversity. The Biomphalaria snails in China have a high genetic similarity with B. straminea sampled from eastern coastal areas of Brazil, which may have originated from the eastern coastal areas of Brazil.
［摘要］ 目的 了解我国藁杆双脐螺来源, 为我国曼氏血吸虫病流行风险评估和双脐螺控制提供依据。方法 选择我 国深圳市观澜河, 大沙河, 深圳水库, 葵涌河上、下游, 新圳河作为采样点, 每个采样点采集10个双脐螺样本, 提取螺样本 基因组DNA。自南美洲巴西米纳斯吉拉斯州、帕拉州、联邦区、伯南布哥州、圣保罗州的5个采样点获得15个藁杆双脐 螺DNA样本。对上述 DNA 样本的细胞色素c 氧化酶亚基 I (cytochrome c oxidase I, COI) 和线粒体16S核糖体RNA (16S ribosomal RNA, 16S rRNA) 基因进行扩增和测序。同时, 从GenBank中下载藁杆双脐螺COI 和16S rRNA 基因序列, 并获取 其采样点信息。将所有基因序列进行比对并构建进化树, 分析我国和南美洲双脐螺样本的遗传相似度和谱系进化关系。结果 从我国双脐螺样本中共获得60个长度为529 bp的COI序列, 其中3个为单倍型。从GenBank中获得165条藁杆双 脐螺COI 序列, 与上述60条序列比对后, 共获得33个单倍型。进化树分析显示, 采自我国的双脐螺3个单倍型聚在一支, 其中单倍型China11与GenBank中获得的来自巴西的3个样本属于同一单倍型。地理进化分析结果显示, 巴西东部沿海 3个采样点的样本有与我国 China11 相同的单倍型, 另2个采样点的样本与 China11 亲缘关系较近。扩增我国双脐螺样本 16S rDNA 基因, 共获得60条长度约为322 bp的序列和2个单倍型。从GenBank中获得70条藁杆双脐螺 16S rDNA 序列。 系统进化树分析显示, 我国双脐螺样本聚为一支, 其中单倍型China64与来自巴西的 229BS 为同一单倍型。将GenBank 中获取的来自巴西南部 25 个采样点的49个藁杆双脐螺 16S rDNA 序列纳入分析, 发现其中3个采样点的藁杆双脐螺与我 国China64有相同的单倍型。综合分析藁杆双脐螺 COI 和 16S rRNA 的地理系统进化关系, 发现仅巴西东部沿海地带样本 与我国双脐螺样本在两个基因片段序列上具有相同单倍型。结论 我国双脐螺为藁杆双脐螺, 遗传多样性较低, 与来 自巴西东部沿海地区的样本遗传学相似度很高, 可能来源于巴西东部沿海地区。.

BACKGROUND: TBX6, a member of the T-box gene family, encodes the transcription factor box 6 that is critical for somite segmentation in vertebrates. It is known that the compound heterozygosity of disruptive variants in trans with a common hypomorphic risk haplotype (T-C-A) in the TBX6 gene contribute to 10% of congenital scoliosis (CS) cases. The deletion of chromosome 17q12 is a rare cytogenetic abnormality, which often leads to renal cysts and diabetes mellitus. However, the affected individuals often exhibit clinical heterogeneity and incomplete penetrance.
METHODS: We here present a Chinese fetus who was shown to have CS by ultrasound examination at 17 weeks of gestation. Trio whole-exome sequencing (WES) was performed to investigate the underlying genetic defects of the fetus. In vitro functional experiments, including western-blotting and luciferase transactivation assay, were performed to determine the pathogenicity of the novel variant of TBX6.
RESULTS: WES revealed the fetus harbored a compound heterozygous variant of c.338_340del (p.Ile113del) and the common hypomorphic risk haplotype of the TBX6 gene. In vitro functional study showed the p.Ile113del variant had no impact on TBX6 expression, but almost led to complete loss of its transcriptional activity. In addition, we identified a 1.85 Mb deletion on 17q12 region in the fetus and the mother. Though there is currently no clinical phenotype associated with this copy number variation in the fetus, it can explain multiple renal cysts in the pregnant woman.
CONCLUSIONS: This study is the first to report a Chinese fetus with a single amino acid deletion variant and a T-C-A haplotype of TBX6. The clinical heterogeneity of 17q12 microdeletion poses significant challenges for prenatal genetic counseling. Our results once again suggest the complexity of prenatal genetic diagnosis.

Mule duck is vitally important to the production of global duck meat. Here, we present two high-quality haplotypes of a female mule duck (haplotype 1 (H1):1.28 Gb, haplotype 2 (H2): 1.40 Gb). The continuity (H1: contig N50 = 14.90 Mb, H2: contig N50 = 15.70 Mb) and completeness (BUSCO: H1 = 96.9%, H2 = 97.3%) are substantially better than those of other duck genomes. We detected the structural variations (SVs) in H1 and H2. We observed a positive correlation between autosome length and the number of SVs. Z chromosome was some deficient in deletions and insertions, but W chromosome was some excessive. A total of 1,451 genes were haplotype specific expression (HSEs). Among them, 737 specifically expressed in H1, and 714 specifically expressed in H2. We found that H1 and H2 HSEs tended to be involved in similar biological processes, such as myometrial relaxation and contraction pathways, muscle structure development and phosphorylation. Our haplotype-resolved genome assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in mule duck.

OBJECTIVE: To analyze a Chinese pedigree with a recombination occurring between the HLA-A/C loci in both parents.
METHODS: A patient who was planning to undergo hematopoietic stem cell transplantation due to \"aplastic anemia\" in February 2022 was selected as the study subject. Peripheral blood samples were collected from the patient, his parents and brother. HLA-A/C/B/DRB1/DQB1 high-resolution typing was carried out by using sequence-based typing and sequence-specific oligonucleotides. The recombination was identified by pedigree analysis. The HLA haplotype of each individual was identified by genealogical analysis. The parentage possibility was determined by short tandem repeat analysis. HLA-A/C/B/DRB1/DRB345/DQA1/DQB1/DPA1/DPB1 were determined with next-generation high-throughput sequence-based typing. The recombination sites were analyzed by family study.
RESULTS: The high parentage possibilities of the family was confirmed by short tandem repeat analysis. Recombination was found between the HLA-A*24:02 A*33:03/C*14:03 in the paternally transmitted haplotype, whilst HLA-A*01:01 A*03:01/C*08:02 was found in the maternally transmitted haplotype, which had resulted in two novel HLA haplotypes in the proband.
CONCLUSIONS: A rare case with simultaneous recombination of the paternal and maternal HLA-A/C loci has been discovered, which may facilitate further study of the mechanisms of the HLA recombination.

Faba bean is an important legume crop. The genetic diversity among faba bean genotypes is very important for the genetic improvement of target traits. A set of 128 fab bean genotypes that are originally from Egypt were used in this study to investigate the genetic diversity and population structure. The 128 genotypes were genotyped using the Single Primer Enrichment Technology (SPET) by which a set of 6759 SNP markers were generated after filtration. The SNP markers were distributed on all chromosomes with a range extending from 822 (Chr. 6) to 1872 (Chr.1). The SNP markers had wide ranges of polymorphic information content (PIC), gene diversity (GD), and minor allele frequency. The analysis of population structure divided the Egyptian faba bean population into five subpopulations. Considerable genetic distance was found among all genotypes, ranging from 0.1 to 0.4. The highly divergent genotype was highlighted in this study and the genetic distance among genotypes ranged from 0.1 and 0.6. Moreover, the structure of linkage disequilibrium was studied, and the analysis revealed a low level of LD in the Egyptian faba bean population. A slow LD decay at the genomic and chromosomal levels was observed. Interestingly, the distribution of haplotype blocks was presented in each chromosome and the number of haplotype block ranged from 65 (Chr. 4) to 156 (Chr. 1). Migration and genetic drift are the main reasons for the low LD in the Egyptian faba bean population. The results of this study shed light on the possibility of the genetic improvement of faba bean crop in Egypt and conducting genetic association analyses to identify candidate genes associated with target traits (e.g. protein content, grain yield, etc.) in this panel.

BACKGROUND: Variations in the innate and adaptive immune response systems are linked to variations in the severity of COVID-19. Natural killer cell (NK) function is regulated by sophisticated receptor system including Killer-cell immunoglobulin-like receptor (KIR) family. We aimed to investigate the impact of possessing certain KIR genes and genotypes on COVID19 severity in Iranians. KIR genotyping was performed on 394 age/sex matched Iranians with no underlying conditions who developed mild and severe COVID- 19. The presence and/or absence of 11 KIR genes were determined using the PCR with sequence specific primers (PCR-SSP).
RESULTS: Patients with mild symptoms had higher frequency ofKIR2DS1 (p = 0.004) and KIR2DS2 (p = 0.017) genes compared to those with severe disease. While KIR3DL3 and deleted variant of KIR2DS4 occurred more frequently in patients who developed a severe form of the disease. In this study, a significant increase of and B haplotype was observed in the Mild group compared to the Severe group (respectively, p = 0.002 and p = 0.02). Also, the prevalence of haplotype A was significantly higher in the Severe group than in the Mild group (p = 0.02).
CONCLUSIONS: These results suggest that the KIR2DS1, KIR2DS, and B haplotype maybe have a protective effect against COVID-19 severity. The results also suggest the inhibitory gene KIR2DL3 and haplotype A are risk factors for the severity of COVID-19.