Annotation of 1350 Common Genetic Variants of the 19 ALDH Multigene Family from Global Human Genome Aggregation Database (gnomAD)

Human aldehyde dehydrogenase (ALDH) is a multigene family with 19 functional members encoding a class of diverse but important enzymes for detoxification or biotransformation of different endogenous and exogenous aldehyde substrates. Genetic mutations in the ALDH genes can cause the accumulation of toxic aldehydes and abnormal carbonyl metabolism and serious human pathologies. However, the physiological functions and substrate specificity of many ALDH genes are still unknown. Although many genetic variants of the ALDH gene family exist in human populations, their phenotype or clinical consequences have not been determined. Using the most comprehensive global human Genome Aggregation Database, gnomAD, we annotated here 1350 common variants in the 19 ALDH genes. These 1350 common variants represent all known genetic polymorphisms with a variant allele frequency of ≥0.1% (or an expected occurrence of ≥1 carrier per 500 individuals) in any of the seven major ethnic groups recorded by gnomAD. We detailed 13 types of DNA sequence variants, their genomic positions, SNP ID numbers, and allele frequencies among the seven major ethnic groups worldwide for each of the 19 ALDH genes. For the 313 missense variants identified in the gnomAD, we used two software algorithms, Polymorphism Phenotyping (PolyPhen) and Sorting Intolerant From Tolerant (SIFT), to predict the consequences of the variants on the structure and function of the enzyme. Finally, gene constraint analysis was used to predict how well genetic mutations were tolerated by selection forces for each of the ALDH genes in humans. Based on the ratio of observed and expected variant numbers in gnomAD, the three ALDH1A gene members, ALDH1A1, ALDH1A2, and ALDH1A3, appeared to have the lowest tolerance for loss-of-function mutations as compared to the other ALDH genes (# observed/# expected ratio 0.15–0.26). These analyses suggest that the ALDH1A1, ALDH1A2, and ALDH1A3 enzymes may serve a more essential function as compared with the other ALDH enzymes; functional loss mutations are much less common in healthy human populations than expected. This informatic analysis may assist the research community in determining the physiological function of ALDH isozymes and associate common variants with clinical phenotypes.

Evolutionary Transition From Water to Land in Vertebrates Illuminated by Basal Ray-Finned Fish Vomeronasal Type 2 Receptor (OlfC) Genes

The vomeronasal type 2 receptor (V2R, also called OlfC) multigene family is found in a broad range of jawed vertebrates from cartilaginous fish to tetrapods. V2Rs encode receptors for food-related amino acids in teleost fish, whereas for peptide pheromones in mammals. In addition, V2Rs of teleost fish are phylogenetically distinct from those of tetrapods, implying a drastic change in the V2Rrepertoire during terrestrial adaptation. To understand the process of diversification of V2Rs in vertebrates from “fish-type” to “tetrapod-type”, we conducted an exhaustive search for V2Rs in cartilaginous fish (chimeras, sharks, and skates) and basal ray-finned fish (reedfish, sterlet, and spotted gar), and compared them with those of teleost, coelacanth, and tetrapods. Phylogenetic and synteny analyses on 1897V2Rs revealed that basal ray-finned fish possess unexpectedly higher number of V2Rs compared with cartilaginous fish, implying that V2Rgene repertoires expanded in the common ancestor of Osteichthyes. Furthermore, reedfish and sterlet possessed various V2Rs that belonged to both “fish-type” and “tetrapod-type”, suggesting that the common ancestor of Osteichthyes possess “tetrapod-type” V2Rs although they inhabited underwater environments. Thus, the unexpected diversity of V2Rs in basal ray-finned fish illuminates the process of how the osteichthyan ancestors adapt from water to land.

The birth-and-death evolution of cytochrome P450 genes in bees

The birth-and-death model of multigene family evolution describes how families can expand by duplication and contract by gene deletion and formation of pseudogenes. The phylogenetic stability of a gene is thought to be related to the degree of functional importance. However, it is unclear how much evolution of a gene in a gene family is driven by adaptive versus neutral processes. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Bees have a high toxin exposure due to their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some bees. Here, we describe the P450s of the orchid bee Euglossa dilemma. Orchid bees are a neotropical clade in which males form perfume bouquets used in courtship displays by collecting a diverse set of volatile compounds, resulting in high chemical compound exposure. We conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We do not find a relationship between the ecology of a bee species and its P450 repertoire. Our analyses reveal that P450 clades can be classified into stable and unstable clades, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures, with signals of adaptive evolution detected, suggesting that both gene duplication and positive selection driving sequence divergence have played a role in the diversification of bee P450s. Our works highlights the complexity of multigene family evolution which does not always follow generalized predictions.

Fishing for Genes: How the Largest Gene Family in the Mammalian Genome was Found (and Why Idiosyncrasy in Exploration Matters)

In 1991, Linda Buck and Richard Axel identified the multigene family expressing odor receptors. Their discovery transformed research on olfaction overnight, and Buck and Axel were awarded the 2004 Nobel Prize in Physiology or Medicine. Behind this success lies another, less visible study about the methodological ingenuity of Buck. This hidden tale holds the key to answering a fundamental question in discovery analysis: What makes specific discovery tools fit their tasks? Why do some strategies turn out to be more fruitful than others? The fit of a method with an experimental system often establishes the success of a discovery. However, the underlying reasoning of discovery is hard to codify. These difficulties point toward an element of discovery analysis routinely sidelined as a mere biographical element in the philosophical analysis of science: the individual discoverer’s role. I argue that the individual researcher is not a replaceable epistemic element in discovery analysis. This article draws on contemporary oral history, including interviews with Buck and other actors key to developments in late 1980s olfaction.

Repeat-Driven Generation of Antigenic Diversity in a Major Human Pathogen, Trypanosoma cruzi

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure.

Classification and characterization of multigene family proteins of African swine fever viruses

African swine fever virus (ASFV) poses serious threats to the pig industry. The multigene family (MGF) proteins are extensively distributed in ASFVs and are generally classified into five families, including MGF-100, MGF-110, MGF-300, MGF-360 and MGF-505. Most MGF proteins, however, have not been well characterized and classified within each family. To bridge this gap, this study first classified MGF proteins into 31 groups based on protein sequence homology and network clustering. A web server for classifying MGF proteins was established and kept available for free at http://www.computationalbiology.cn/MGF/home.html. Results showed that MGF groups of the same family were most similar to each other and had conserved sequence motifs; the genetic diversity of MGF groups varied widely, mainly due to the occurrence of indels. In addition, the MGF proteins were predicted to have large structural and functional diversity, and MGF proteins of the same MGF family tended to have similar structure, location and function. Reconstruction of the ancestral states of MGF groups along the ASFV phylogeny showed that most MGF groups experienced either the copy number variations or the gain-or-loss changes, and most of these changes happened within strains of the same genotype. It is found that the copy number decrease and the loss of MGF groups were much larger than the copy number increase and the gain of MGF groups, respectively, suggesting the ASFV tended to lose MGF proteins in the evolution. Overall, the work provides a detailed classification for MGF proteins and would facilitate further research on MGF proteins.

Genome-wide identification, and phylogenetic and expression profiling analyses, of XTH gene families in Brassica rapa L. and Brassica oleracea L.

Background Xyloglucan endotransglucosylase/hydrolase genes (XTHs) are a multigene family and play key roles in regulating cell wall extensibility in plant growth and development. Brassica rapa and Brassica oleracea contain XTHs, but detailed identification and characterization of the XTH family in these species, and analysis of their tissue expression profiles, have not previously been carried out. Results In this study, 53 and 38 XTH genes were identified in B. rapa and B. oleracea respectively, which contained some novel members not observed in previous studies. All XTHs of B. rapa, B. oleracea and Arabidopsis thaliana could be classified into three groups, Group I/II, III and the Early diverging group, based on phylogenetic relationships. Gene structures and motif patterns were similar within each group. All XTHs in this study contained two characteristic conserved domains (Glyco_hydro and XET_C). XTHs are located mainly in the cell wall but some are also located in the cytoplasm. Analyses of the mechanisms of gene family expansion revealed that whole-genome triplication (WGT) events and tandem duplication (TD) may have been the major mechanisms accounting for the expansion of the XTH gene family. Interestingly, TD genes all belonged to Group I/II, suggesting that TD was the main reason for the largest number of genes being in these groups. B. oleracea had lost more of the XTH genes, the conserved domain XET_C and the conserved active-site motif EXDXE compared with B. rapa, consistent with asymmetrical evolution between the two Brassica genomes. A majority of XTH genes exhibited different tissue-specific expression patterns based on RNA-seq data analyses. Moreover, there was differential expression of duplicated XTH genes in the two species, indicating that their functional differentiation occurred after B. rapa and B. oleracea diverged from a common ancestor. Conclusions We carried out the first systematic analysis of XTH gene families in B. rapa and B. oleracea. The results of this investigation can be used for reference in further studies on the functions of XTH genes and the evolution of this multigene family.