Genomic Regions Research Articles

The biogeographic and evolutionary processes shaping population divergence in Laupala.

Speciation generates biodiversity and the mechanisms involved are thought to vary across the tree of life and across environments. For example, well-studied adaptive radiations are thought to be fuelled by divergent ecological selection, but additionally are influenced heavily by biogeographic, genomic and demographic factors. Mechanisms of non-adaptive radiations, producing ecologically cryptic taxa, have been less well-studied but should likewise be influenced by these latter factors. Comparing among contexts can help pinpoint universal mechanisms and outcomes, especially if we integrate biogeographic, ecological and evolutionary processes. We investigate population divergence in the swordtail cricket Laupala cerasina, a wide-spread endemic on Hawai'i Island and one of 38 ecologically cryptic Laupala species. The nine sampled populations show striking population genetic structure at small spatio-temporal scales. The rapid differentiation among populations and species of Laupala shows that neither a specific geographical context nor ecological opportunity are pre-requisites for rapid divergence. Spatio-temporal patterns in population divergence, population size change, and gene flow are aligned with the chronosequence of the four volcanoes on which L. cerasina occurs and reveal the composite effects of geological dynamics and Quaternary climate change on population dynamics. Spatio-temporal patterns in genetic variation along the genome reveal the interplay of genetic and genomic architecture in shaping population divergence. In early phases of divergence, we find elevated differentiation in genomic regions harbouring mating song loci. In later stages of divergence, we find a signature of linked selection that interacts with recombination rate variation. Comparing our findings with recent work on complementary systems supports the conclusion that mostly universal factors influence the speciation process.

Metabolic Profile of the Genome-Reduced Bacillus subtilis Strain IIG-Bs-27-39: An Attractive Chassis for Recombinant Protein Production.

The Gram-positive bacterium Bacillus subtilis is extensively used in the industry for the secretory production of proteins with commercial value. To further improve its performance, this microbe has been the subject of extensive genome engineering efforts, especially the removal of large genomic regions that are dispensable or even counterproductive. Here, we present the genome-reduced B. subtilis strain IIG-Bs-27-39, which was obtained through systematic deletion of mobile genetic elements, as well as genes for extracellular proteases, sporulation, flagella formation, and antibiotic production. Different from previously characterized genome-reduced B. subtilis strains, the IIG-Bs-27-39 strain was still able to grow on minimal media. We used this feature to benchmark strain IIG-Bs-27-39 against its parental strain 168 with respect to heterologous protein production and metabolic parameters during bioreactor cultivation. The IIG-Bs-27-39 strain presented superior secretion of difficult-to-produce staphylococcal antigens, as well as higher specific growth rates and biomass yields. At the metabolic level, changes in byproduct formation and internal amino acid pools were observed, whereas energetic parameters such as the ATP yield, ATP/ADP levels, and adenylate energy charge were comparable between the two strains. Intriguingly, we observed a significant increase in the total cellular NADPH level during all tested conditions and increases in the NAD+ and NADP(H) pools during protein production. This indicates that the IIG-Bs-27-39 strain has more energy available for anabolic processes and protein production, thereby providing a link between strain physiology and production performance. On this basis, we conclude that the genome-reduced strain IIG-Bs-27-39 represents an attractive chassis for future biotechnological applications.

Genetic Control of Tolerance to Drought Stress in Wild Soybean (Glycine soja) at the Vegetative and the Germination Stages

Drought stress, which is becoming more prevalent due to climate change, is a significant abiotic factor that adversely impacts crop production and yield stability. Cultivated soybean (Glycine max), a versatile crop for humans and animals, exhibits sensitivity to drought, resulting in reduced growth and development under drought conditions. However, few genetic studies have assessed wild soybean’s (Glycine soja) response to drought stress. In this work, we conducted a genome-wide association study (GWAS) and analysis of wild soybean accessions to identify loci responsible for drought tolerance at the vegetative (n = 187) and the germination stages (n = 135) using the available resequencing data. The GWAS analysis of the leaf wilting score (LWS) identified eight single-nucleotide polymorphisms (SNPs) on chromosomes 10, 11, and 19. Of these, wild soybeans with both SNPs on chromosomes 10 (adenine) and 11 (thymine) produced lower LWS, indicating that these SNPs have an important role in the genetic effect on LWS for drought tolerance at the vegetative stage. At the germination stage, nine SNPs associated with five phenotypic measurements were identified on chromosomes 6, 9, 10, 13, 16, and 17, and the genomic regions identified at the germination stage were different from those identified for the LWS, supporting our previous finding that there may not be a robust correlation between the genes influencing phenotypes at the germination and vegetative stages. This research will benefit marker-assisted breeding programs aimed at enhancing drought tolerance in soybeans.

Evolution at Spike protein position 519 in SARS-CoV-2 facilitated adaptation to humans

As the COVID-19 pandemic enters its fourth year, the pursuit of identifying a progenitor virus to SARS-CoV-2 and understanding the mechanism of its emergence persists, albeit against the backdrop of intensified efforts to monitor the ongoing evolution of the virus and the influx of new mutations. Surprisingly, few residues hypothesized to be essential for SARS-CoV-2 emergence and adaptation to humans have been validated experimentally, despite the importance that these mutations could contribute to the development of effective antivirals. To remedy this, we searched for genomic regions in the SARS-CoV-2 genome that show evidence of past selection around residues unique to SARS-CoV-2 compared with closely related coronaviruses. In doing so, we identified a residue at position 519 in Spike within the receptor binding domain that holds a static histidine in human-derived SARS-CoV-2 sequences but an asparagine in SARS-related coronaviruses from bats and pangolins. In experimental validation, the SARS-CoV-2 Spike protein mutant carrying the putatively ancestral H519N substitution showed reduced replication in human lung cells, suggesting that the histidine residue contributes to viral fitness in the human host. Structural analyses revealed a potential role of Spike residue 519 in mediating conformational transitions necessary for Spike prior to binding with ACE2. Pseudotyped viruses bearing the putatively ancestral N519 also demonstrated significantly reduced infectivity in cells expressing the human ACE2 receptor compared to H519. ELISA data corroborated that H519 enhances Spike binding affinity to the human ACE2 receptor compared to the putatively ancestral N519. Collectively, these findings suggest that the evolutionary transition at position 519 of the Spike protein played a critical role in SARS-CoV-2 emergence and adaptation to the human host. Additionally, this residue presents as a potential drug target for designing small molecule inhibitors tailored to this site.

Genome-Wide Methylation Landscape Uncovers the Role of DNA Methylation in Ramie (Boehmeria nivea L.) Bast Fiber Growth

Ramie is one of the most important fiber crops in China, with fibers extracted from stem barks having been used as textile materials for thousands of years. DNA methylation is an important epigenetic modification involved in plant growth and development. However, the role of methylation in ramie fiber growth remains poorly understood. In the present study, we investigated the DNA methylation landscape of the nuclear genome in bark sections taken from the top (TPS) and the middle (MPS) of the stems of ramie plants, which represent different stages of fiber growth, using whole-genome bisulfite sequencing. We detected 7,709,555 and 8,508,326 5-methylcytosines in the TPS and MPS genomes, respectively. The distribution of methylation across three sequence contexts, CG, CHG, and CHH, varied greatly among gene elements, with methylation at CHH being the most prevalent. Comparison of methylation levels between the TPS and MPS genomes revealed 23.162 Mb of differentially methylated genomic regions, encompassing 9485 genes. Among these differentially methylated genes, 841 exhibited altered expression in the MPS genome. Notably, an SND2 ortholog Bni05G006779 showed a negative correlation between its expression and methylation levels. Overexpression of Bni05G006779 in Arabidopsis dramatically increased the number of xylem fibers and the secondary wall thickness of the fibers in the stems of transgenic plants. These findings provide important insights into the involvement of DNA methylation in regulating ramie fiber growth.

Genome-wide association study reveals marker-trait associations for major agronomic traits in proso millet (Panicum miliaceum L.).

The pilot-scale genome-wide association study in the US proso millet identified twenty marker-trait associations for five morpho-agronomic traits identifying genomic regions for future studies (e.g. molecular breeding and map-based cloning). Proso millet (Panicum miliaceum L.) is an ancient grain recognized for its excellent water-use efficiency and short growing season. It is an indispensable part of the winterwheat-based dryland cropping system in the High Plains of the USA. Its grains are endowed with high nutritional and health-promoting properties, making it increasingly popular in the global market for healthy grains. There is a dearth of genomic resources in proso millet for developing molecular tools to complement conventional breeding for developing high-yielding varieties. Genome-wide association study (GWAS) is a widely used method to dissect the genetics of complex traits. In this pilot study of the first-ever GWAS in the US proso millet, 71 globally diverse genotypes of 109 the US proso millet core collection were evaluated for five major morpho-agronomic traits at two locations in western Nebraska, and GWAS was conducted to identify single nucleotide polymorphisms (SNPs) associated with these traits. Analysis of variance showed that there was a significant difference among the genotypes, and all five traits were also found to be highly correlated with each other. Sequence reads from genotyping-by-sequencing (GBS) were used to identify 11,147 high-quality bi-allelic SNPs. Population structure analysis with those SNPs showed stratification within the core collection. The GWAS identified twenty marker-trait associations (MTAs) for the five traits. Twenty-nine putative candidate genes associated with the five traits were also identified. These genomic regions can be used to develop genetic markers for marker-assisted selection in proso millet breeding.

Congenital syndromic Chiari-like malformation (CSCM) in Holstein cattle: towards unravelling of possible genetic causes

BackgroundChiari malformation type II (CMII) was originally reported in humans as a rare disorder characterized by the downward herniation of the hindbrain and towering cerebellum. The congenital brain malformation is usually accompanied by spina bifida, a congenital spinal anomaly resulting from incomplete closure of the dorsal aspect of the spinal neural tube, and occasionally by other lesions. A similar disorder has been reported in several animal species, including cattle, particularly as a congenital syndrome. A cause of congenital syndromic Chiari-like malformation (CSCM) in cattle has not been reported to date. We collected a series of 14 CSCM-affected Holstein calves (13 purebred, one Red Danish Dairy F1 cross) and performed whole-genome sequencing (WGS). WGS was performed on 33 cattle, including eight cases with parents (trio-based; group 1), three cases with one parent (group 2), and three single cases (solo-based; group 3).ResultsSequencing-based genome-wide association study of the 13 Holstein calves with CSCM and 166 controls revealed no significantly associated genome region. Assuming a single Holstein breed-specific recessive allele, no region of shared homozygosity was detected suggesting heterogeneity. Subsequent filtering for protein-changing variants that were only homozygous in the genomes of the individual cases allowed the identification of two missense variants affecting different genes, SHC4 in case 4 in group 1 and WDR45B in case 13 in group 3. Furthermore, these two variants were only observed in Holstein cattle when querying WGS data of > 5,100 animals. Alternatively, potential de novo mutational events were assessed in each case. Filtering for heterozygous private protein-changing variants identified one DYNC1H1 frameshift variant as a candidate causal dominant acting allele in case 12 in group 3. Finally, the presence of larger structural DNA variants and chromosomal abnormalities was investigated in all cases. Depth of coverage analysis revealed two different partial monosomies of chromosome 2 segments in cases 1 and 7 in group 1 and a trisomy of chromosome 12 in the WDR45B homozygous case 13 in group 3.ConclusionsThis study presents for the first time a detailed genomic evaluation of CSCM in Holstein cattle and suggests an unexpected genetic and allelic heterogeneity considering the mode of inheritance, as well as the type of variant. For the first time, we propose candidate causal variants that may explain bovine CSCM in a certain proportion of affected calves. We present cattle as a large animal model for human CMII and propose new genes and genomic variants as possible causes for related diseases in both animals and humans.

Genomic insights into the shared and distinct genetic architecture of cognitive function and schizophrenia

Cognitive impairment is a major determinant of functional outcomes in schizophrenia, however, understanding of the biological mechanisms underpinning cognitive dysfunction in the disorder remains incomplete. Here, we apply Genomic Structural Equation Modelling to identify latent cognitive factors capturing genetic liabilities to 12 cognitive traits measured in the UK Biobank. We identified three broad factors that underly the genetic correlations between the cognitive tests. We explore the overlap between latent cognitive factors, schizophrenia, and schizophrenia symptom dimensions using a complementary set of statistical approaches, applied to data from the latest schizophrenia genome-wide association study (Ncase = 53,386, Ncontrol = 77,258) and the Thematically Organised Psychosis study (Ncase = 306, Ncontrol = 1060). Global genetic correlations showed a significant moderate negative genetic correlation between each cognitive factor and schizophrenia. Local genetic correlations implicated unique genomic regions underlying the overlap between schizophrenia and each cognitive factor. We found substantial polygenic overlap between each cognitive factor and schizophrenia and biological annotation of the shared loci implicated gene-sets related to neurodevelopment and neuronal function. Lastly, we show that the common genetic determinants of the latent cognitive factors are not predictive of schizophrenia symptoms in the Norwegian Thematically Organized Psychosis cohort. Overall, these findings inform our understanding of cognitive function in schizophrenia by demonstrating important differences in the shared genetic architecture of schizophrenia and cognitive abilities.

Structural Basis for the Interaction Between Yeast Chromatin Assembly Factor 1 and Proliferating Cell Nuclear Antigen

Proliferating cell nuclear antigen (PCNA), the homotrimeric eukaryotic sliding clamp protein, recruits and coordinates the activities of a multitude of proteins that function on DNA at the replication fork. Chromatin assembly factor 1 (CAF-1), one such protein, is a histone chaperone that deposits histone proteins onto DNA immediately following replication. The interaction between CAF-1 and PCNA is essential for proper nucleosome assembly at silenced genomic regions. Most proteins that bind PCNA contain a PCNA-interacting peptide (PIP) motif, a conserved motif containing only eight amino acids. Precisely how PCNA is able to discriminate between binding partners at the replication fork using only these small motifs remains unclear. Yeast CAF-1 contains a PIP motif on its largest subunit, Cac1. We solved the crystal structure of the PIP motif of CAF-1 bound to PCNA using a new strategy to produce stoichiometric quantities of one PIP motif bound to each monomer of PCNA. The PIP motif of CAF-1 binds to the hydrophobic pocket on the front face of PCNA in a similar manner to most known PIP-PCNA interactions. However, several amino acids immediately flanking either side of the PIP motif bind the IDCL or C-terminus of PCNA, as observed for only a couple other known PIP-PCNA interactions. Furthermore, mutational analysis suggests positively charged amino acids in these flanking regions are responsible for the low micromolar affinity of CAF-1 for PCNA, whereas the presence of a negative charge upstream of the PIP prevents a more robust interaction with PCNA. These results provide additional evidence that positive charges within PIP-flanking regions of PCNA-interacting proteins are crucial for specificity and affinity of their recruitment to PCNA at the replication fork.

Epigenomic features associated with body temperature stabilize tissues during cold exposure in cold-resistant pigs

Cold stress in low-temperature environments can trigger changes in gene expression, but epigenomics regulation of temperature stability in vital tissues, including the fat and diencephalon, is still unclear. Here, we explore the cold-induced changes in epigenomic features in the diencephalon and fat tissues of two cold-resistant Chinese pig breeds, Min and Enshi black (ES) pigs, utilizing H3K27ac CUT&Tag, RNA-seq, and selective signature analysis. Our results show significant alterations in H3K27ac modifications in the diencephalon of Min pigs and the fat of ES pigs after cold exposure. Dramatic changes in H3K27ac modifications in Min pigs are primarily associated with genes involved in energy metabolism and hormone regulation, whereas those in ES pigs are primarily associated with immunity-related genes. Moreover, transcription factors PRDM1 and HSF1, which show evidence of selection, are enriched in genomic regions presenting cold-responsive alterations in H3K27ac modification in the Min pig diencephalon and ES pig fat, respectively. Our results indicate the diversity of epigenomic response mechanisms to cold exposure between Min and ES pigs, providing unique epigenetic resources for studies of low-temperature adaptation in large mammals.

Genetic analyses point to alterations in immune-related pathways underpinning the association between psychiatric disorders and COVID-19.

Current literature suggests that people with psychiatric disorders have a higher risk of COVID-19 infection and a worse prognosis of the disease. We aimed to study the genetic contribution to these associations across seven psychiatric disorders as well as a general psychopathology factor (P-factor) and determine whether these are unique or shared across psychiatric disorders using statistical genetic techniques. Using the largest available genome-wide association studies (GWAS), we found a significant genetic overlap between depression, ADHD, PTSD, and the P-factor with both COVID-19 infection and hospitalization, and between anxiety and COVID-19 hospitalization. We used pairwise GWAS to examine this overlap on a fine-grained scale and identified specific regions of the genome shared between several psychiatric disorders, the P-factor, and COVID-19. Gene-based analysis in these genomic regions suggested possible links with immune-related pathways such as thyroid homeostasis, inflammation, and stress response. Finally, we show preliminary evidence for causal associations between depression, ADHD, PTSD, and the P-factor, and higher COVID-19 infection and hospitalization using Mendelian Randomization and Latent Causal Variable methods. Our results support the hypothesis that the relationship between psychiatric disorders and COVID-19 risk is likely due to shared alterations in immune-related pathways and is not a result of environmental factors alone, shedding light on potentially viable therapeutic targets.

The Identification of the Peanut Wild Relative Arachis stenosperma as a Source of Resistance to Stem Rot and Analyses of Genomic Regions Conferring Disease Resistance through QTL Mapping

Peanut stem rot, also known as white mold, poses a significant threat to peanut production. It is typically managed using fungicides and moderately resistant cultivars. Cultivars with higher resistance can reduce fungicide dependency and increase sustainability. This study explores the potential of wild peanut species in stem rot resistance breeding programs by enhancing genetic diversity in cultivated peanut. Through greenhouse and field evaluations, 13 allotetraploid hybrids with Arachis stenosperma as one of the parents showed superior resistance compared to other wild genotypes. The genomic regions that confer the stem rot resistance were further identified by genotyping and phenotyping an F2 population derived from the allotetraploid ValSten1 (A. valida × A. stenosperma)4× and A. hypogaea cv. TifGP-2. A linkage map was constructed from 1926 SNP markers. QTL analysis revealed both beneficial and deleterious loci, with two resistance-associated QTLs derived from A. stenosperma and four susceptibility loci, two from A. stenosperma and two from A. valida. This is the first study that evaluated peanut-compatible wild-derived allotetraploids for stem rot resistance and that identified wild-derived QTLs for resistance to this pathogen. The allotetraploid hybrid ValSten1, that has A. stenosperma as one of the parents, offers a resource for resistance breeding. Markers associated with resistance QTLs can facilitate introgression from ValSten1 into cultivated peanut varieties in future breeding efforts, potentially reducing reliance on chemical control measures.

Identifying genomic regions associated with C4 photosynthetic activity and leaf anatomy in Alloteropsis semialata.

C4 photosynthesis is a complex trait requiring multiple developmental and metabolic alterations. Despite this complexity, it has independently evolved over 60 times. However, our understanding of the transition to C4 is complicated by the fact that variation in photosynthetic type is usually segregated between species that diverged a long time ago. Here, we perform a genome-wide association study (GWAS) using the grass Alloteropsis semialata, the only known species to have C3, intermediate, and C4 accessions that recently diverged. We aimed to identify genomic regions associated with the strength of the C4 cycle (measured using δ13C), and the development of C4 leaf anatomy. Genomic regions correlated with δ13C include regulators of C4 decarboxylation enzymes (RIPK), nonphotochemical quenching (SOQ1), and the development of Kranz anatomy (SCARECROW-LIKE). Regions associated with the development of C4 leaf anatomy in the intermediate individuals contain additional leaf anatomy regulators, including those responsible for vein patterning (GSL8) and meristem determinacy (GIF1). The parallel recruitment of paralogous leaf anatomy regulators between A. semialata and other C4 lineages implies the co-option of these genes is context-dependent, which likely has implications for the engineering of the C4 trait into C3 species.

Comprehensive survey of conserved RNA secondary structures in full-genome alignment of Hepatitis C virus

Hepatitis C virus (HCV) is a plus-stranded RNA virus that often chronically infects liver hepatocytes and causes liver cirrhosis and cancer. These viruses replicate their genomes employing error-prone replicases. Thereby, they routinely generate a large ‘cloud’ of RNA genomes (quasispecies) which—by trial and error—comprehensively explore the sequence space available for functional RNA genomes that maintain the ability for efficient replication and immune escape. In this context, it is important to identify which RNA secondary structures in the sequence space of the HCV genome are conserved, likely due to functional requirements. Here, we provide the first genome-wide multiple sequence alignment (MSA) with the prediction of RNA secondary structures throughout all representative full-length HCV genomes. We selected 57 representative genomes by clustering all complete HCV genomes from the BV-BRC database based on k-mer distributions and dimension reduction and adding RefSeq sequences. We include annotations of previously recognized features for easy comparison to other studies. Our results indicate that mainly the core coding region, the C-terminal NS5A region, and the NS5B region contain secondary structure elements that are conserved beyond coding sequence requirements, indicating functionality on the RNA level. In contrast, the genome regions in between contain less highly conserved structures. The results provide a complete description of all conserved RNA secondary structures and make clear that functionally important RNA secondary structures are present in certain HCV genome regions but are largely absent from other regions. Full-genome alignments of all branches of Hepacivirus C are provided in the supplement.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

IntroductionEffective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.MethodsWe investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.ResultsSnap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes

BackgroundThe 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.ResultsWe 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).ConclusionsIn 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.

The association of rs17713054 with Neanderthal origin at 3p21.31 locus with the severity of COVID-19 in Iranian patients

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 (Pvalue = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (Pvalue = 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.

Emergence of the ZNF675 Gene During Primate Evolution-Influenced Human Neurodevelopment Through Changing HES1 Autoregulation.

In this study, we investigated recurrent copy number variations (CNVs) in the 19p12 locus, which are associated with neurodevelopmental disorders. The two genes in this locus, ZNF675 and ZNF681, arose via gene duplication in primates, and their presence in several pathological CNVs in the human population suggests that either or both of these genes are required for normal human brain development. ZNF675 and ZNF681 are members of the Krüppel-associated box zinc finger (KZNF) protein family, a class of transcriptional repressors important for epigenetic silencing of specific genomic regions. About 170 primate-specific KZNFs are present in the human genome. Although KZNFs are primarily associated with repressing retrotransposon-derived DNA, evidence is emerging that they can be co-opted for other gene regulatory processes. We show that genetic deletion of ZNF675 causes developmental defects in cortical organoids, and our data suggest that part of the observed neurodevelopmental phenotype is mediated by a gene regulatory role of ZNF675 on the promoter of the neurodevelopmental gene Hes family BHLH transcription factor 1 (HES1). We also find evidence for the recently evolved regulation of genes involved in neurological disorders, microcephalin 1 and sestrin 3. We show that ZNF675 interferes with HES1 auto-inhibition, a process essential for the maintenance of neural progenitors. As a striking example of how some KZNFs have integrated into preexisting gene expression networks, these findings suggest the emergence of ZNF675 has caused a change in the balance of HES1 autoregulation. The association of ZNF675 CNV with human developmental disorders and ZNF675-mediated regulation of neurodevelopmental genes suggests that it evolved into an important factor for human brain development.

First identification of canine parvovirus -2a/2b variant in unvaccinated domestic dogs with gastrointestinal signs in Türkiye.

Canine parvovirus type 2 (CPV-2) is the most common enteric virus that infects canids. CPV is the causative agent of a contagious disease defined mostly by clinical gastrointestinal signs in dogs. During the late 1970s, CPV-2 emerged as a new virus capable of infecting domestic dogs and growing across the world. The VP2 gene stands out as a key determinant in the pathogenicity, antigenicity, and host interactions of CPV-2. The molecular characterization of the VP2 gene is crucial for understanding CPV evolution and epidemiology. Genes encoding the VP2 protein were sequenced and compared to reference strains worldwide. The maximum likelihood method was used to build a phylogenetic tree using CPV VP2 gene nucleotide sequences. Our phylogenetic analysis of the VP2 gene revealed that five strains were very similar and clustered together, and three strains were in the 2b clade, whereas the other two were in the 2a/2b clade. This paper reports the molecular characterization of two novel CPV-2a/2b subtypes in dogs with gastrointestinal symptoms. Genetic analysis was conducted on a CPV genomic region encompassing one of the open reading frames (ORFs) encoding the structural protein VP2. Sequence analysis indicates new and unreported sequence changes, mainly affecting the VP2 gene, which includes the mutations Ser297Ala and Leu87Met. This study represents the first evidence of a new CPV-2a/2b subtype in Türkiye. Due to VP2's crucial role in encoding the capsid protein of CPV-2 and its significant involvement in the host-virus interaction, it is critical to closely monitor its evolutionary changes and be cautious while searching for novel or pre-existing subtypes. This study highlights the significance of continuous molecular research for acquiring more insights on the circulation of novel CPV mutants.

Prediction and validation of putative candidate genes underlying drought‐tolerant QTLs through in silico and RT‐PCR approaches in rice (Oryza sativa)

AbstractIn the present study, the aim was to predict and validate putative candidate genes underlying drought‐tolerant quantitative trait loci (QTLs) in rice crop using in silico approaches and real‐time polymerase chain reaction (RT‐PCR). The genes underlying major drought‐tolerant QTLs which have been reported by data mining, sequence variation, gene ontology analysis, quantitative traits gene finder and gene expression analysis were subjected to RiceVarmap software to design primers, and only a few variants gave the SNP/InDel primers; thus, finally, 15 primers were ultimately selected, which were used in identification of differentially expressed genes (DEGs) among contrasting rice genotypes IR 64 and N 22 for drought tolerance trait using quantitative RT‐PCR studies by providing drought stress treatment during panicle initiation stage. In this investigation, we predicted 11 genes as candidate genes underlying drought‐tolerant QTLs. Out of these, only four QTLs were found responsible for the major effect in drought tolerance regions such as QTL‐Qsn‐4b, QTL‐rn7a, QTL‐Qtgw‐2a and QTL‐phc4.1 and 11 prioritized candidates were identified that expressed in leaf tissues. Only four primers belong to two QTLs, primer vg0712623096 from QTL‐rn7a (LOC_Os07g22450) located on chromosome‐7 encoding NAC domain‐containing protein and the primers vg0431750843(LOC_Os04g53310) encoding soluble starch synthase 3‐ chloroplast precursor, vg0432626757 (LOC_Os04g54850) encoding pectin acetylesterase domain‐containing protein and vg0433031562 (LOC_Os04g55520) encoding AP2 domain‐containing protein, from QTL‐Qsn‐4b, located on chromosome‐4 found to have higher differential expression in N 22 in comparison with IR 64 during drought stress as per quantitative RT‐PCR 2–ΔΔCt values. Considering the overall study, these four primers/genes were identified as candidate genes underlying genomic regions governing drought tolerance. Therefore, these putative candidate genes could be focussed for further functional analysis to exploit in rice breeding.