Effects of Population Bottleneck and Balancing Selection on the Chinese Alligator Are Revealed by Locus-Specific Characterization of MHC Genes

Balancing selection shapes population differentiation of major histocompatibility complex genes in wild golden snub-nosed monkeys

Small and isolated populations face several intrinsic risks, such as genetic drift, inbreeding depression, and reduced gene flow. Thus, patterns of genetic diversity and differentiation have become an important focus of conservation genetics research. The golden snub-nosed monkey Rhinopithecus roxellana, an endangered species endemic to China, has experienced rapid reduction in population size and severe population fragmentation over the past few decades. We measured the patterns of genetic diversity and population differentiation using both neutral microsatellites and adaptive major histocompatibility complex (MHC) genes in 2 R. roxellana populations (DPY and GNG) distributed on the northern and southern slopes of the Qinling Mountains, respectively. Eight MHC-linked haplotypes formed by 5 DQA1 alleles, 5 DQB1 alleles, 5 DRB1 alleles, and 4 DRB2 alleles were detected in the 2 populations. The larger GNG population showed higher genetic variation for both MHC and microsatellites than the smaller DPY population, suggesting an effect of genetic drift on genetic variation. Genetic differentiation index (F ST) outlier analyses, principal coordinate analysis (PCoA), and inferred population genetic structure showed lower genetic differentiation in the MHC variations than microsatellites, suggesting that pathogen-mediated balancing selection, rather than local adaptation, homogenized the MHC genes of both populations. This study indicates that both balancing selection and genetic drift may shape genetic variation and differentiation in small and fragmented populations.

Copy number and selection of MHC genes in ruminants are related to habitat, average life span and diet

The ruminants, as the main group of livestock, have been extensively studied in terms of their physiology, endocrinology, biochemistry, genetics, and nutrition. Despite the wide geographic distribution and habitat diversity of animals in this group, their ecology and evolution remain poorly understood. In this study, we analyzed the gene copy number, selection, and ecological and evolutionary processes that have affected the evolution of major histocompatibility complex (MHC) genes across ruminant lineages based on available genomic data. The 51 species analyzed represented all six families of ruminants. Our finding indicated that the architecture of the MHC region is conserved in ruminants, but with variable copy numbers of MHC-I, MHC-IIA, and MHC-IIB genes. No lineage-specific gene duplication was observed in the MHC genes. The phylogenetic generalized least squares regression (PGLS) model revealed association between ecological and biological factors (habitat and lifespan) and gene duplication in DQA and DQB, but not in DRB. The selection pressure of DQA and DQB were related with lifespan, diet, and the ratio of genetic repeat elements. These results suggest that the MHC evolution in ruminants, including copy number and selection, has been influenced by genetic repeat elements, pathogen exposure risk, and intrinsic cost of possessing multiple MHC genes.

Harnessing γ-TMT Genetic Variations and Haplotypes for Vitamin E Diversity in the Korean Rice Collection

Gamma-tocopherol methyltransferase (γ-TMT), a key gene in the vitamin E biosynthesis pathway, significantly influences the accumulation of tocochromanols, thereby determining rice nutritional quality. In our study, we analyzed the γ-TMT gene in 475 Korean rice accessions, uncovering 177 genetic variants, including 138 SNPs and 39 InDels. Notably, two functional SNPs, tmt-E2-28,895,665-G/A and tmt-E4-28,896,689-A/G, were identified, causing substitutions from valine to isoleucine and arginine to glycine, respectively, across 93 accessions. A positive Tajima's D value in the indica group suggests a signature of balancing selection. Haplotype analysis revealed 27 haplotypes, with two shared between cultivated and wild accessions, seven specific to cultivated accessions, and 18 unique to wild types. Further, profiling of vitamin E isomers in 240 accessions and their association with haplotypes revealed that Hap_2, distinguished by an SNP in the 3' UTR (tmt-3UTR-28,897,360-T/A) exhibited significantly lower α-tocopherol (AT), α-tocotrienol (AT3), total tocopherol, and total tocotrienol, but higher γ-tocopherol (GT) in the japonica group. Additionally, in the indica group, Hap_2 showed significantly higher AT, AT3, and total tocopherol, along with lower GT and γ-tocotrienol, compared to Hap_19, Hap_20, and Hap_21. Overall, this study highlights the genetic landscape of γ-TMT and provides a valuable genetic resource for haplotype-based breeding programs aimed at enhancing nutritional profiles.

Global conservation prioritization areas in three dimensions of crocodilian diversity

Crocodilians are a taxonomic group of large predators with important ecological and evolutionary benefits for ecosystem functioning in the face of global change. Anthropogenic actions affect negatively crocodilians' survival and more than half of the species are threatened with extinction worldwide. Here, we map and explore three dimensions of crocodilian diversity on a global scale. To highlight the ecological importance of crocodilians, we correlate the spatial distribution of species with the ecosystem services of nutrient retention in the world. We calculate the effectiveness of global protected networks in safeguarding crocodilian species and provide three prioritization models for conservation planning. Our results show the main hotspots of ecological and evolutionary values are in southern North, Central and South America, west-central Africa, northeastern India, and southeastern Asia. African species have the highest correlation to nutrient retention patterns. Twenty-five percent of the world's crocodilian species are not significantly represented in the existing protected area networks. The most alarming cases are reported in northeastern India, eastern China, and west-central Africa, which include threatened species with low or non-significant representation in the protected area networks. Our highest conservation prioritization model targets southern North America, east-central Central America, northern South America, west-central Africa, northeastern India, eastern China, southern Laos, Cambodia, and some points in southeastern Asia. Our research provides a global prioritization scheme to protect multiple dimensions of crocodilian diversity for achieving effective conservation outcomes.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Major histocompatibility complex genomic investigation of endangered Chinese alligator provides insights into the evolution of tetrapod major histocompatibility complex and survival of critically bottlenecked species

Background

The major histocompatibility complex (MHC) gene family, a vital immune gene family in vertebrates, helps animals defend against pathogens. The polymorphism of MHC genes is important for a species and is considered to be caused by the numerous alleles of MHC antigen-presenting genes. However, the mechanism of this process is unclear due to the lack of data on the MHC structure. The evolutionary trajectories of the tetrapod MHC are also unclear because of insufficient studies on the reptile MHC architecture. Here, we studied the Chinese alligator (Alligator sinensis), which experienced a population bottleneck, but the population increased rapidly in the past 30 years and is proposed to have a unique MHC system to face pathogenic challenges.

Results

We successfully constructed a 2 Mb MHC region using bacterial artificial chromosome (BAC) library and genome data of the Chinese alligator and checked the antigen-presenting genes using transcriptome data and the rapid amplification of cDNA ends (RACE) technique. The MHC architecture reported here uncovers adjacent Class I and Class II subregions and a unique CD1 subregion. This newly added information suggested that the Class I-II structure pattern was more ancient in tetrapods and helped reconstruct the evolution of the MHC region architecture. We also found multiple groups of MHC class I (MHC-I) (12 duplicated loci, belonging to three groups, two of which were novel) and MHC class II (MHC-II) (11 duplicated loci, belonging to two groups) inside the 2 Mb MHC region, and there were three more duplicated MHC-I loci outside it. These highly duplicated antigen-presenting loci had differences in expression, amino acid length of antigen-presenting exons, and splice signal of exon and intron, which together promoted the polymorphism of duplicated genes. Although the MHC antigen-presenting genes were identified as monomorphic or oligomorphic in our previous population study, the loci with high copy numbers and many differences can make up for this loss, presenting another mechanism for polymorphism in antigen presentation. These MHC-I and MHC-IIB loci with low polymorphism for each locus, but high numbers in all, may also contribute to MHC antigen-presenting binding variability in a population.

Conclusion

To summarize, the fine MHC region architecture of reptiles presented in this study completes the evolutionary trajectories of the MHC structure in tetrapods, and these distinctive MHC gene groups in the Chinese alligator may have helped this species to expand rapidly in the past recent years.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Genetic diversity of major histocompatibility complex class I genes in Zootoca vivipara

The Major Histocompatibility Complex (MHC), as a family of highly polymorphic genes associated with immunity in the genome of the vertebrate, has become an important indicator for assessing the evolutionary potential of wildlife. In order to better protect Zootoca vivipara in the Greater Khingan Range and Lesser Khingan Range, to understand the genetic structure of Z. vivipara, and to explore the mechanism and phylogenetic relationship of the gene polymorphisms, the MHC molecular marker method was used to analyze Z. vivipara population. Forty-seven alleles were obtained from four populations. The four populations were highly polymorphic, rich in genetic information, and had significant genetic diversity. There were certain inbreeding phenomena. There was a high degree of genetic differentiation among populations, which was caused by genetic drift and natural selection. The sequence undergoes genetic duplication and recombination. The existence of trans-species polymorphism was found in the constructed phylogenetic tree. The present study provides a theoretical basis for species protection of Z. vivipara.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.