Evolutionary paths to macrolide resistance in a Neisseria commensal converge on ribosomal genes through short sequence duplications

Neisseria commensals are an indisputable source of resistance for their pathogenic relatives. However, the evolutionary paths commensal species take to reduced susceptibility in this genus have been relatively underexplored. Here, we leverage in vitro selection as a powerful screen to identify the genetic adaptations that produce azithromycin resistance (≥ 2 μg/mL) in the Neisseria commensal, N. elongata. Across multiple lineages (n = 7/16), we find mutations that reduce susceptibility to azithromycin converge on the locus encoding the 50S ribosomal L34 protein (rpmH) and the intergenic region proximal to the 30S ribosomal S3 protein (rpsC) through short tandem duplication events. Interestingly, one of the laboratory evolved mutations in rpmH is identical (7LKRTYQ12), and two nearly identical, to those recently reported to contribute to high-level azithromycin resistance in N. gonorrhoeae. Transformations into the ancestral N. elongata lineage confirmed the causality of both rpmH and rpsC mutations. Though most lineages inheriting duplications suffered in vitro fitness costs, one variant showed no growth defect, suggesting the possibility that it may be sustained in natural populations. Ultimately, studies like this will be critical for predicting commensal alleles that could rapidly disseminate into pathogen populations via allelic exchange across recombinogenic microbial genera.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

Whole-Genome Duplication and Host Genotype Affect Rhizosphere Microbial Communities

Plants influence the composition of their associated microbial communities, yet the underlying host-associated genetic determinants are typically unknown. Genome duplication events are common in the evolutionary history of plants and affect many plant traits.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Genome-Wide Prediction, Functional Divergence, and Characterization of Stress-Responsive BZR Transcription Factors in B. napus

BRASSINAZOLE RESISTANT (BZR) are transcriptional factors that bind to the DNA of targeted genes to regulate several plant growth and physiological processes in response to abiotic and biotic stresses. However, information on such genes in Brassica napus is minimal. Furthermore, the new reference Brassica napus genome offers an excellent opportunity to systematically characterize this gene family in B. napus. In our study, 21 BnaBZR genes were distributed across 19 chromosomes of B. napus and clustered into four subgroups based on Arabidopsis thaliana orthologs. Functional divergence analysis among these groups evident the shifting of evolutionary rate after the duplication events. In terms of structural analysis, the BnaBZR genes within each subgroup are highly conserved but are distinctive within groups. Organ-specific expression analyses of BnaBZR genes using RNA-seq data and quantitative real-time polymerase chain reaction (qRT-PCR) revealed complex expression patterns in plant tissues during stress conditions. In which genes belonging to subgroups III and IV were identified to play central roles in plant tolerance to salt, drought, and Sclerotinia sclerotiorum stress. The insights from this study enrich our understanding of the B. napus BZR gene family and lay a foundation for future research in improving rape seed environmental adaptability.

Genome-Wide Characterization and Abiotic Stresses Expression Analysis of Annexin Family Genes in Poplar

Poplar is an illustrious industrial woody plant with rapid growth, providing a range of materials, and having simple post-treatment. Various kinds of environmental stresses limit its output. Plant annexin (ANN) is a calcium-dependent phospholipid-binding protein involved in plant metabolism, growth and development, and cooperatively regulating drought resistance, salt tolerance, and various stress responses. However, the features of the PtANN gene family and different stress responses remain unknown in poplar. This study identified 12 PtANN genes in the P. trichocarpa whole-genome and PtANNs divided into three subfamilies based on the phylogenetic tree. The PtANNs clustered into the same clade shared similar gene structures and conserved motifs. The 12 PtANN genes were located in ten chromosomes, and segmental duplication events were illustrated as the main duplication method. Additionally, the PtANN4 homogenous with AtANN1 was detected localized in the cytoplasm and plasma membrane. In addition, expression levels of PtANNs were induced by multiple abiotic stresses, which indicated that PtANNs could widely participate in response to abiotic stress. These results revealed the molecular evolution of PtANNs and their profiles in response to abiotic stress.

A chromosome-scale genome assembly for the holly (Ilex polyneura) provides insights into genomic adaptations to elevation in Southwest China

Southwest China is a plant diversity hotspot. The near-cosmopolitan genus Ilex (c. 664 spp., Aquifoliaceae) reaches its maximum diversity in this region, with many narrow-range and a few widespread species. Divergent selection on widespread species leads to local adaptation, with consequences for both conservation and utilization, but is counteracted by geneflow. Many Ilex species are utilized as teas, medicines, ornamentals, honey plants, and timber, but variation below the species level is largely uninvestigated. We therefore studied the widespread Ilex polyneura, which occupies most of the elevational range available and is cultivated for its decorative leafless branches with persistent red fruits. We assembled a chromosome-scale genome using approximately 100x whole genome long-read and short-read sequencing combined with Hi-C sequencing. The genome is approximately 727.1 Mb, with a contig N50 size of 5 124 369 bp and a scaffold N50 size of 36 593 620 bp, for which the BUSCO score was 97.6%, and 98.9% of the assembly was anchored to 20 pseudochromosomes. Out of 32 838 genes predicted, 96.9% were assigned functions. Two whole genome duplication events were identified. Using this genome as a reference, we conducted a population genomics study of 112 individuals from 21 populations across the elevation range using restriction site-associated DNA sequencing (RADseq). Most populations clustered into four clades separated by distance and elevation. Selective sweep analyses identified 34 candidate genes potentially under selection at different elevations, with functions related to responses to abiotic and biotic stresses. This first high-quality genome in the Aquifoliales will facilitate the further domestication of the genus.

Genome-wide Analysis of CNVs in Three Populations of Tibetan Sheep using Whole-genome Resequencing

Background: Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptations in the Qinghai-Tibetan Plateau.Results: A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected in three Tibetan sheep populations. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis at each CNVR showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between HTS and VTS, and RUNX1 and LOC101111988 were significantly differentiated between VTS and OTS. Meaningfully, the duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in QTP, which deserves further research in detail.Conclusions: In this study, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

Ankfn1-mutant vestibular defects require loss of both ancestral and derived paralogs for penetrance in zebrafish

How and to what degree gene duplication events create regulatory innovation, redundancy, or neofunctionalization remain important questions in animal evolution and comparative genetics. Ankfn1 genes are single copy in most invertebrates, partially duplicated in jawed vertebrates, and only the derived copy retained in most mammals. Null mutations in the single mouse homolog have vestibular and neurological abnormalities. Null mutation of the single Drosophila homolog is typically lethal with severe sensorimotor deficits in rare survivors. The functions and potential redundancy of paralogs in species with two copies is not known. Here we define a vestibular role for Ankfn1 homologs in zebrafish based on simultaneous disruption of each locus. Zebrafish with both paralogs disrupted showed vestibular defects and early lethality from swim bladder inflation failure. One intact copy at either locus was sufficient to prevent major phenotypes. Our results show that vertebrate Ankfn1 genes are required for vestibular-related functions, with at least partial redundancy between ancestral and derived paralogs.

Molecular evidence for segmental duplication across chromosomes of soybean using transcription factor gene family

Duplication of genome is an important genetic innovation. Large genome size (1.1 Gb) along with ancient and recent duplication events make the soybean genome more complex. Analyzing the distribution and duplication event in soybean transcription family genes, the segmental duplication within chromosomes was revealed. Our study provides a strong evidence that the large segmental duplication event in genome architecture and evolution of soybean genome using simple method of sequence and order analysis of TF genes. Finally, a scheme for interrelationship of different chromosomes has been proposed.