Chromosome-level genome assembly and demographic history of Euryodendron excelsum in monotypic genus endemic to China

Euryodendron excelsum is in a monotypic genus Euryodendron, endemic to China. It has intermediate morphisms in the Pentaphylacaceae or Theaceae families, which make it distinct. Due to anthropogenic disturbance, E. excelsum is currently found in very restricted and fragmented areas with extremely small populations. Although much research and effort has been applied towards its conservation, its long-term survival mechanisms and evolutionary history remain elusive, especially from a genomic aspect. Therefore, using a combination of long/short whole genome sequencing, RNA sequencing reads, and Hi-C data, we assembled and annotated a high-quality genome for E. excelsum. The genome assembly of E. excelsum comprised 1,059,895,887 bp with 99.66% anchored into 23 pseudo-chromosomes and a 99.0% BUSCO completeness. Comparative genomic analysis revealed the expansion of terpenoid and flavonoid secondary metabolite genes, and displayed a tandem and/or proximal duplication framework of these genes. E. excelsum also displayed genes associated with growth, development, and defence adaptation from whole genome duplication. Demographic analysis indicated that its fluctuations in population size and its recent population decline were related to cold climate changes. The E. excelsum genome assembly provides a highly valuable resource for evolutionary and ecological research in the future, aiding its conservation, management, and restoration.

Origin and early divergence of tandem duplicated sorbitol transporter genes in Rosaceae: insights from evolutionary analysis of the SOT gene family in angiosperms

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol transporters (SOTs) play a vital role in facilitating sorbitol allocation from source to sink organs and sugar accumulation in sink organs. While prior research has addressed gene duplications within the SOT gene family in Rosaceae, the precise origin and evolutionary dynamics of these duplications remain unclear, largely due to the complicated interplay of whole genome duplications and tandem duplications. Here, we investigated the synteny relationships among all identified Polyol/Monosaccharide Transporter (PLT) genes in 61 angiosperm genomes and SOT genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of PLTs and SOTs across diverse plant lineages. We found that Rosaceae SOTs, as PLT family members, originated from a pair of tandemly duplicated PLT genes within Class III-A. Furthermore, our investigation highlights the role of lineage-specific and synergistic duplications in Amygdaloideae in contributing to the expansion of SOTs in Rosaceae plants. Collectively, our findings provide insights into the genomic origins, duplication events, and subsequent divergence of SOT gene family members. Such insights lay a crucial foundation for comprehensive functional characterizations in future studies.

Genome-wide identification and evolutionary analysis of the FGF gene family in buffalo

Fibroblast growth factors (FGFs) are important polypeptide growth factors that play a critical role in many developmental processes, including differentiation, cell proliferation, and migration in mammals. This study employs in silico analyses to characterize the FGF gene family in buffalo, investigating their genome-wide identification, physicochemical properties, and evolutionary patterns. For this purpose, genomic and proteomic sequences of buffalo, cattle, goat, and sheep were retrieved from NCBI database. We identified a total of 22 FGF genes in buffalo. Physicochemical properties observed through ProtParam tool showed notable features of these proteins including in-vitro instability, thermostability, hydrophilicity, and basic nature. Phylogenetic analysis grouped 22 identified genes into nine sub-families based on evolutionary relationships. Additionally, analysis of gene structure, motif patterns, and conserved domains using TBtools revealed the remarkable conservation of this gene family across selected species throughout the course of evolution. Comparative amino acid analysis performed through ClustalW demonstrated significant conservation between buffalo and cattle FGF proteins. Mutational analysis showed three non-synonymous mutations at positions R103 > G, P7 > L, and E98 > Q in FGF4, FGF6, and FGF19, respectively in buffalo. Duplication events revealed only one segmental duplication (FGF10/FGF22) in buffalo and two in cattle (FGF10/FGF22 and FGF13/FGF13-like) with Ka/Ks values <1 indicating purifying selection pressure for these duplications. Comparison of protein structures of buffalo, goat, and sheep exhibited more similarities in respective structures. In conclusion, our study highlights the conservation of the FGF gene family in buffalo during evolution. Furthermore, the identified non-synonymous mutations may have implications for the selection of animals with better performance.Communicated by Ramaswamy H. Sarma.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

The pathogenic p.Gln319Ter variant is not causing congenital adrenal hyperplasia when inherited in one of the duplicated CYP21A2 genes

Objective

The study aimed to identify the pathogenic status of p.Gln319Ter (NM_000500.7: c.955C>T) variant when inherited in a single CYP21A2 gene (bimodular RCCX haplotype) and to discriminate between a non-causing congenital adrenal hyperplasia (CAH) allele when inherited in a duplicated and functional CYP21A2 gene context (trimodular RCCX haplotype).

Methods

38 females and 8 males with hyperandrogenemia, previously screened by sequencing and identified as carriers for the pathogenic p.Gln319Ter, were herein tested by multiplex ligation-dependent probe amplification (MLPA) and a real-time PCR Copy number Variation (CNV) assay.

Results

Both MLPA and real-time PCR CNV analyses confirmed a bimodular and pathogenic RCCX haplotype with a single CYP21A2 in 19/46 (41.30%) p.Gln319Ter carriers and who in parallel all shared elevated 17-OHP levels. The remaining 27 individuals that also carried the p.Gln319Ter exhibited low 17-OHP levels as a result of their carriership of a duplicated CYP21A2 with a trimodular RCCX haplotype. Interestingly, all of these individuals also carried in linkage disequilibrium with p.Gln319Ter two single nucleotide polymorphisms, the c.293-79G>A (rs114414746) in intron 2 and the c.*12C>T (rs150697472) in the 3'-UTR. Therefore, these variants can be used to distinguish between pathogenic and non-pathogenic genomic contexts of the c.955T (p.Gln319) in the genetic diagnosis of congenital adrenal hyperplasia (CAH).

Conclusion

The employed methodologies identified a considerable number of individuals with non-pathogenic p.Gln319Ter from the individuals that typically carry the pathogenic p.Gln319Ter in a single CYP21A2. Therefore, it is extremely important the detection of such haplotypes for the prenatal diagnosis, treatment and genetic counseling in patients with CAH.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Genome-wide analysis of R2R3-MYB genes in cultivated peanut (Arachis hypogaea L.): Gene duplications, functional conservation, and diversification

The cultivated Peanut (Arachis hypogaea L.), an important oilseed and edible legume, are widely grown worldwide. The R2R3-MYB transcription factor, one of the largest gene families in plants, is involved in various plant developmental processes and responds to multiple stresses. In this study we identified 196 typical R2R3-MYB genes in the genome of cultivated peanut. Comparative phylogenetic analysis with Arabidopsis divided them into 48 subgroups. The motif composition and gene structure independently supported the subgroup delineation. Collinearity analysis indicated polyploidization, tandem, and segmental duplication were the main driver of the R2R3-MYB gene amplification in peanut. Homologous gene pairs between the two subgroups showed tissue specific biased expression. In addition, a total of 90 R2R3-MYB genes showed significant differential expression levels in response to waterlogging stress. Furthermore, we identified an SNP located in the third exon region of AdMYB03-18 (AhMYB033) by association analysis, and the three haplotypes of the SNP were significantly correlated with total branch number (TBN), pod length (PL) and root-shoot ratio (RS ratio), respectively, revealing the potential function of AdMYB03-18 (AhMYB033) in improving peanut yield. Together, these studies provide evidence for functional diversity in the R2R3-MYB genes and will contribute to understanding the function of R2R3-MYB genes in peanut.

Differential Expression of Two Copies of the irmA Gene in the Enteroaggregative E. coli Strain 042

Gene duplications significantly impact the gene repertoires of both eukaryotic and prokaryotic microorganisms. The genomes of pathogenic Escherichia coli strains share a group of duplicated genes whose function is mostly unknown. The irmA gene is one of the duplicates encoded in several pathogenic E. coli strains. The function of its gene product was investigated in the uropathogenic E. coli strain CFT073, which contains a single functional copy. The IrmA protein structure mimics that of human interleukin receptors and likely plays a role during infection. The enteroaggregative E. coli strain 042 contains two functional copies of the irmA gene. In the present work, we investigated their biological roles. The irmA_4509 allele is expressed under several growth conditions. Its expression is modulated by the global regulators OxyR and Hha, with optimal expression at 37°C and under nutritional stress conditions. Expression of the irmA_2244 allele can only be detected when the irmA_4509 allele is knocked out. Differences in the promoter regions of both alleles account for their differential expression. Our results show that under several environmental conditions, the expression of the IrmA protein in strain 042 is dictated by the irmA_4509 allele. The irmA_2244 allele appears to play a backup role to ensure IrmA expression when the irmA_4509 allele loses its function. IMPORTANCE Gene duplications occur in prokaryotic genomes at a detectable frequency. In many instances, the biological function of the duplicates is unknown, and hence, the significance of the presence of multiple copies of these genes remains unclear. In pathogenic E. coli isolates, the irmA gene can be present either as a single copy or in two or more copies. We focused our work on studying why a different pathogenic E. coli strain encodes two functional copies of the irmA gene. We show that under several environmental conditions, one of the alleles dictates IrmA expression, and the second remains silent. The latter allele is only expressed when the former is silenced. The presence of more than one functional copy of the irmA gene in some pathogenic E. coli strains can result in sufficient expression of this virulence factor during the infection process.

Gene duplication to the Y chromosome in Trinidadian Guppies

Differences in allele frequencies at autosomal genes between males and females in a population can result from two scenarios. First, unresolved sexual conflict over survival can produce allelic differentiation between the sexes. However, given the substantial mortality costs required to produce allelic differences between males and females at each generation, it remains unclear how many loci within the genome experience significant sexual conflict over survival. Alternatively, recent studies have shown that similarity between autosomal and Y sequences can create perceived allelic differences between the sexes. However, Y duplications are most likely in species with large non-recombining regions, in part because they simply represent larger targets for duplications. We assessed the genomes of 120 wild-caught guppies, which experience extensive predation- and pathogen-induced mortality and have a relatively small ancestral Y chromosome. We identified seven autosomal genes that show allelic differences between male and female adults. Five of these genes show clear evidence of whole or partial gene duplication between the Y chromosome and the autosomes. The remaining two genes show evidence of partial homology to the Y. Overall, our findings suggest that the guppy genome experiences a very low level of unresolved sexual conflict over survival, and instead the Y chromosome, despite its small ancestral size and recent origin, may nonetheless accumulate genes with male-specific functions.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Plant Metabolic Gene Clusters: Evolution, Organization, and Their Applications in Synthetic Biology

Plants are a remarkable source of high-value specialized metabolites having significant physiological and ecological functions. Genes responsible for synthesizing specialized metabolites are often clustered together for a coordinated expression, which is commonly observed in bacteria and filamentous fungi. Similar to prokaryotic gene clustering, plants do have gene clusters encoding enzymes involved in the biosynthesis of specialized metabolites. More than 20 gene clusters involved in the biosynthesis of diverse metabolites have been identified across the plant kingdom. Recent studies demonstrate that gene clusters are evolved through gene duplications and neofunctionalization of primary metabolic pathway genes. Often, these clusters are tightly regulated at nucleosome level. The prevalence of gene clusters related to specialized metabolites offers an attractive possibility of an untapped source of highly useful biomolecules. Accordingly, the identification and functional characterization of novel biosynthetic pathways in plants need to be worked out. In this review, we summarize insights into the evolution of gene clusters and discuss the organization and importance of specific gene clusters in the biosynthesis of specialized metabolites. Regulatory mechanisms which operate in some of the important gene clusters have also been briefly described. Finally, we highlight the importance of gene clusters to develop future metabolic engineering or synthetic biology strategies for the heterologous production of novel metabolites.

Gene duplication at the Fascicled ear1 locus controls the fate of inflorescence meristem cells in maize

The maize ear is unbranched and terminates in a single point. The ear and tassel inflorescences of Fascicled ear mutants fail to grow as a single point and instead are branched. This phenotype results from the misexpression of duplicated transcription factors, ZMM8 and DRL2. We hypothesize that these gene rearrangements create regulatory sequences that cause misexpression in early inflorescence meristems, thus activating a laminar program, ablating the meristem, and producing branches. This work demonstrates that zmm8 and drl2 must be restricted from the inflorescence meristem to maintain its terminal point, and conversely, a mechanism by which branching may be imposed. Manipulation of these genes can be used to alter plant architecture, potentially to improve agronomic traits.

Plant meristems are self-renewing groups of pluripotent stem cells that produce lateral organs in a stereotypical pattern. Of interest is how the radially symmetrical meristem produces laminar lateral organs. Both the male and female inflorescence meristems of the dominant Fascicled ear (Fas1) mutant fail to grow as a single point and instead show deep branching. Positional cloning of two independent Fas1 alleles identified an ∼160 kb region containing two floral genes, the MADS-box gene, zmm8, and the YABBY gene, drooping leaf2 (drl2). Both genes are duplicated within the Fas1 locus and spatiotemporally misexpressed in the mutant inflorescence meristems. Increased zmm8 expression alone does not affect inflorescence development; however, combined misexpression of zmm8, drl2, and their syntenic paralogs zmm14 and drl1, perturbs meristem organization. We hypothesize that misexpression of the floral genes in the inflorescence and their potential interaction cause ectopic activation of a laminar program, thereby disrupting signaling necessary for maintenance of radially symmetrical inflorescence meristems. Consistent with this hypothesis, RNA sequencing and in situ analysis reveal altered expression patterns of genes that define distinct zones of the meristem and developing leaf. Our findings highlight the importance of strict spatiotemporal patterns of expression for both zmm8 and drl2 and provide an example of phenotypes arising from tandem gene duplications.

Duplications and Functional Convergence of Intestinal Carbohydrate-Digesting Enzymes

Vertebrate diets and digestive physiologies vary tremendously. Although the contribution of ecological and behavioral features to such diversity is well documented, the roles and identities of individual intestinal enzymes shaping digestive traits remain largely unexplored. Here, we show that the sucrase-isomaltase (SI)/maltase-glucoamylase (MGAM) dual enzyme system long assumed to be the conserved disaccharide and starch digestion framework in all vertebrates is absent in many lineages. Our analyses indicate that independent duplications of an ancestral SI gave rise to the mammalian-specific MGAM, as well as to other duplicates in fish and birds. Strikingly, the duplicated avian enzyme exhibits similar activities to MGAM, revealing an unexpected case of functional convergence. Our results highlight digestive enzyme variation as a key uncharacterized component of dietary diversity in vertebrates.

Evolution and Expression of the Immune System of a Facultatively Anadromous Salmonid

Vertebrates have evolved a complex immune system required for the identification of and coordinated response to harmful pathogens. Migratory species spend periods of their life-cycle in more than one environment, and their immune system consequently faces a greater diversity of pathogens residing in different environments. In facultatively anadromous salmonids, individuals may spend parts of their life-cycle in freshwater and marine environments. For species such as the brown trout Salmo trutta, sexes differ in their life-histories with females more likely to migrate to sea while males are more likely to stay and complete their life-cycle in their natal river. Salmonids have also undergone a lineage-specific whole genome duplication event, which may provide novel immune innovations but our current understanding of the differences in salmonid immune expression between the sexes is limited. We characterized the brown trout immune gene repertoire, identifying a number of canonical immune genes in non-salmonid teleosts to be duplicated in S. trutta, with genes involved in innate and adaptive immunity. Through genome-wide transcriptional profiling (“RNA-seq”) of male and female livers to investigate sex differences in gene expression amplitude and alternative splicing, we identified immune genes as being generally male-biased in expression. Our study provides important insights into the evolutionary consequences of whole genome duplication events on the salmonid immune gene repertoire and how the sexes differ in constitutive immune expression.

Mitochondrial copper and phosphate transporter specificity was defined early in the evolution of eukaryotes

The mitochondrial carrier family protein SLC25A3 transports both copper and phosphate in mammals, yet in Saccharomyces cerevisiae the transport of these substrates is partitioned across two paralogs: PIC2 and MIR1. To understand the ancestral state of copper and phosphate transport in mitochondria, we explored the evolutionary relationships of PIC2 and MIR1 orthologs across the eukaryotic tree of life. Phylogenetic analyses revealed that PIC2-like and MIR1-like orthologs are present in all major eukaryotic supergroups, indicating an ancient gene duplication created these paralogs. To link this phylogenetic signal to protein function, we used structural modeling and site-directed mutagenesis to identify residues involved in copper and phosphate transport. Based on these analyses, we generated an L175A variant of mouse SLC25A3 that retains the ability to transport copper but not phosphate. This work highlights the utility of using an evolutionary framework to uncover amino acids involved in substrate recognition by mitochondrial carrier family proteins.

Evolution of MicroRNA Biogenesis Genes in the Sterlet (Acipenser ruthenus) and Other Polyploid Vertebrates

MicroRNAs play a crucial role in eukaryotic gene regulation. For a long time, only little was known about microRNA-based gene regulatory mechanisms in polyploid animal genomes due to difficulties of polyploid genome assembly. However, in recent years, several polyploid genomes of fish, amphibian, and even invertebrate species have been sequenced and assembled. Here we investigated several key microRNA-associated genes in the recently sequenced sterlet (Acipenser ruthenus) genome, whose lineage has undergone a whole genome duplication around 180 MYA. We show that two paralogs of drosha, dgcr8, xpo1, and xpo5 as well as most ago genes have been retained after the acipenserid-specific whole genome duplication, while ago1 and ago3 genes have lost one paralog. While most diploid vertebrates possess only a single copy of dicer1, we strikingly found four paralogs of this gene in the sterlet genome, derived from a tandem segmental duplication that occurred prior to the last whole genome duplication. ago1,3,4 and exportins1,5 look to be prone to additional segment duplications producing up to four-five paralog copies in ray-finned fishes. We demonstrate for the first time exon microsatellite amplification in the acipenserid drosha2 gene, resulting in a highly variable protein product, which may indicate sub- or neofunctionalization. Paralogous copies of most microRNA metabolism genes exhibit different expression profiles in various tissues and remain functional despite the rediploidization process. Subfunctionalization of microRNA processing gene paralogs may be beneficial for different pathways of microRNA metabolism. Genetic variability of microRNA processing genes may represent a substrate for natural selection, and, by increasing genetic plasticity, could facilitate adaptations to changing environments.

Ancient Evolutionary Origin of Intrinsically Disordered Cancer Risk Regions

Cancer is a heterogeneous genetic disease that alters the proper functioning of proteins involved in key regulatory processes such as cell cycle, DNA repair, survival, or apoptosis. Mutations often accumulate in hot-spots regions, highlighting critical functional modules within these proteins that need to be altered, amplified, or abolished for tumor formation. Recent evidence suggests that these mutational hotspots can correspond not only to globular domains, but also to intrinsically disordered regions (IDRs), which play a significant role in a subset of cancer types. IDRs have distinct functional properties that originate from their inherent flexibility. Generally, they correspond to more recent evolutionary inventions and show larger sequence variations across species. In this work, we analyzed the evolutionary origin of disordered regions that are specifically targeted in cancer. Surprisingly, the majority of these disordered cancer risk regions showed remarkable conservation with ancient evolutionary origin, stemming from the earliest multicellular animals or even beyond. Nevertheless, we encountered several examples where the mutated region emerged at a later stage compared with the origin of the gene family. We also showed the cancer risk regions become quickly fixated after their emergence, but evolution continues to tinker with their genes with novel regulatory elements introduced even at the level of humans. Our concise analysis provides a much clearer picture of the emergence of key regulatory elements in proteins and highlights the importance of taking into account the modular organisation of proteins for the analyses of evolutionary origin.

Young sex chromosomes in plants and animals

A major reason for studying plant sex chromosomes is that they may often be 'young' systems. There is considerable evidence for the independent evolution of separate sexes within plant families or genera, in some cases showing that the maximum possible time during which their sex-determining genes have existed must be much shorter than those of several animal taxa. Consequently, their sex-linked regions could either have evolved soon after genetic sex determination arose or considerably later. Plants, therefore, include species with both young and old systems. I review several questions about the evolution of sex-determining systems and sex chromosomes that require studies of young systems, including: the kinds of mutations involved in the transition to unisexual reproduction from hermaphroditism or monoecy (a form of functional hermaphroditism); the times when they arose; and the extent to which the properties of sex-linked regions of genomes reflect responses to new selective situations created by the presence of a sex-determining locus. I also evaluate which questions are best studied in plants, vs other suitable candidate organisms. Studies of young plant systems can help understand general evolutionary processes that are shared with the sex chromosomes of other organisms.

Genome-Wide Identification and Expression Profiling of the Polygalacturonase (PG) and Pectin Methylesterase (PME) Genes in Grapevine (Vitis vinifera L.)

In pectin regulation, polygalacturonases (PGs) and pectin methylesterases (PMEs) are critical components in the transformation, disassembly network, and remodeling of plant primary cell walls. In the current study, we identified 36 PG and 47 PME genes using the available genomic resources of grapevine. Herein, we provide a comprehensive overview of PGs and PMEs, including phylogenetic and collinearity relationships, motif and gene structure compositions, gene duplications, principal component analysis, and expression profiling during developmental stages. Phylogenetic analysis of PGs and PMEs revealed similar domain composition patterns with Arabidopsis. The collinearity analysis showed high conservation and gene duplications with purifying selection. The type of duplications also varied in terms of gene numbers in PGs (10 dispersed, 1 proximal, 12 tandem, and 13 segmental, respectively) and PMEs (23 dispersed, 1 proximal, 16 tandem, and 7 segmental, respectively). The tissue-specific response of PG and PME genes based on the reported transcriptomic data exhibited diverged expression patterns in various organs during different developmental stages. Among PGs, VvPG8, VvPG10, VvPG13, VvPG17, VvPG18, VvPG19, VvPG20, VvPG22, and VvPG23 showed tissue- or organ-specific expression in majority of the tissues during development. Similarly, in PMEs, VvPME3, VvPME4, VvPME5, VvPME6, VvPME19, VvPME21, VvPME23, VvPME29, VvPME31, and VvPME32 suggested high tissue-specific response. The gene ontology (GO), Kyoto Encyclopedia of Genes and Genomics (KEGG) enrichment, and cis-elements prediction analysis also suggested the putative functions of PGs and PMEs in plant development, such as pectin and carbohydrate metabolism, and stress activities. Moreover, qRT-PCR validation of 32 PG and PME genes revealed their role in various organs of grapevines (i.e., root, stem, tendril, inflorescence, flesh, skins, and leaves). Therefore, these findings will lead to novel insights and encourage cutting-edge research on functional characterization of PGs and PMEs in fruit crop species.

Consequences of CYP2D6 Copy-Number Variation for Pharmacogenomics in Psychiatry

Pharmacogenomics represents a potentially powerful enhancement to the current standard of care for psychiatric patients. However, a variety of biological and technical challenges must be addressed in order to provide adequate clinical decision support for personalized prescribing and dosing based on genomic data. This is particularly true in the case of CYP2D6, a key drug-metabolizing gene, which not only harbors multiple genetic variants known to affect enzyme function but also shows a broad range of copy-number and hybrid alleles in various patient populations. Here, we describe several challenges in the accurate measurement and interpretation of data from the CYP2D6 locus including the clinical consequences of increased copy number. We discuss best practices for overcoming these challenges and then explore various current and future applications of pharmacogenomic analysis of CYP2D6 in psychiatry.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.

Molecular evolution of key metabolic genes during transitions to C4 and CAM photosynthesis

PREMISE OF THE STUDY

Next-generation sequencing facilitates rapid production of well-sampled phylogenies built from very large genetic data sets, which can then be subsequently exploited to examine the molecular evolution of the genes themselves. We present an evolutionary analysis of 83 gene families (19 containing carbon-concentrating mechanism (CCM) genes, 64 containing non-CCM genes) in the portullugo clade (Caryophyllales), a diverse lineage of mostly arid-adapted plants that contains multiple evolutionary origins of all known photosynthesis types in land plants (C₃ , C₄ , CAM, C₄ -CAM, and various intermediates).

METHODS

We inferred a phylogeny of 197 individuals from 167 taxa using coalescent-based approaches and individual gene family trees using maximum likelihood. Positive selection analyses were conducted on individual gene family trees with a mixed effects model of evolution (MEME). We devised new indices to compare levels of convergence and prevalence of particular residues between CCM and non-CCM genes and between species with different photosynthetic pathways.

KEY RESULTS

Contrary to expectations, there were no significant differences in the levels of positive selection detected in CCM versus non-CCM genes. However, we documented a significantly higher level of convergent amino acid substitutions in CCM genes, especially in C₄ taxa.

CONCLUSIONS

Our analyses reveal a new suite of amino acid residues putatively important for C₄ and CAM function. We discuss both the advantages and challenges of using targeted enrichment sequence data for exploratory studies of molecular evolution.

Side-by-side secretion of Late Palaeozoic diverged courtship pheromones in an aquatic salamander

Males of the advanced salamanders (Salamandroidea) attain internal fertilization without a copulatory organ by depositing a spermatophore on the substrate in the environment, which females subsequently take up with their cloaca. The aquatically reproducing modern Eurasian newts (Salamandridae) have taken this to extremes, because most species do not display close physical contact during courtship, but instead largely rely on females following the male track at spermatophore deposition. Although pheromones have been widely assumed to represent an important aspect of male courtship, molecules able to induce the female following behaviour that is the prelude for successful insemination have not yet been identified. Here, we show that uncleaved sodefrin precursor-like factor (SPF) protein pheromones are sufficient to elicit such behaviour in female palmate newts (Lissotriton helveticus). Combined transcriptomic and proteomic evidence shows that males simultaneously tail-fan multiple ca 20 kDa glycosylated SPF proteins during courtship. Notably, molecular dating estimates show that the diversification of these proteins already started in the late Palaeozoic, about 300 million years ago. Our study thus not only extends the use of uncleaved SPF proteins outside terrestrially reproducing plethodontid salamanders, but also reveals one of the oldest vertebrate pheromone systems.

The evolutionary dynamics of major regulators for sexual development among Hymenoptera species

All hymenopteran species, such as bees, wasps and ants, are characterized by the common principle of haplodiploid sex determination in which haploid males arise from unfertilized eggs and females from fertilized eggs. The underlying molecular mechanism has been studied in detail in the western honey bee Apis mellifera, in which the gene complementary sex determiner (csd) acts as primary signal of the sex determining pathway, initiating female development by csd-heterozygotes. Csd arose from gene duplication of the feminizer (fem) gene, a transformer (tra) ortholog, and mediates in conjunction with transformer2 (tra2) sex-specific splicing of fem. Comparative molecular analyses identified fem/tra and its downstream target doublesex (dsx) as conserved unit within the sex determining pathway of holometabolous insects. In this study, we aim to examine evolutionary differences among these key regulators. Our main hypothesis is that sex determining key regulators in Hymenoptera species show signs of coevolution within single phylogenetic lineages. We take advantage of several newly sequenced genomes of bee species to test this hypothesis using bioinformatic approaches. We found evidences that duplications of fem are restricted to certain bee lineages and notable amino acid differences of tra2 between Apis and non-Apis species propose structural changes in Tra2 protein affecting co-regulatory function on target genes. These findings may help to gain deeper insights into the ancestral mode of hymenopteran sex determination and support the common view of the remarkable evolutionary flexibility in this regulatory pathway.

Genomics and genetics of gonadotropin beta-subunit genes: Unique FSHB and duplicated LHB/CGB loci

The follicle stimulating hormone (FSH), luteinizing hormone (LH) and chorionic gonadotropin (HCG) play a critical role in human reproduction. Despite the common evolutionary ancestry and functional relatedness of the gonadotropin hormone beta (GtHB) genes, the single-copy FSHB (at 11p13) and the multi-copy LHB/CGB genes (at 19q13.32) exhibit locus-specific differences regarding their genomic context, evolution, genetic variation and expressional profile. FSHB represents a conservative vertebrate gene with a unique function and it is located in a structurally stable gene-poor region. In contrast, the primate-specific LHB/CGB gene cluster is located in a gene-rich genomic context and demonstrates an example of evolutionary young and unstable genomic region. The gene cluster is shaped by a constant balance between selection that acts on specific functions of the loci and frequent gene conversion events among duplicons. As the transcription of the GtHB genes is rate-limiting in the assembly of respective hormones, the genomic and genetic context of the FSHB and the LHB/CGB genes largely affects the profile of the hormone production.

Tissue-driven hypothesis of genomic evolution and sequence-expression correlations

To maintain normal physiological functions, different tissues may have different developmental constraints on expressed genes. Consequently, the evolutionary tolerance for genomic evolution varies among tissues. Here, we formulate this argument as a "tissue-driven hypothesis" based on the stabilizing selection model. Moreover, several predicted genomic correlations are tested by the human-mouse microarray data. Our results are as follows. First, between the human and mouse, we have elaborated the among-tissue covariation between tissue expression distance (E(ti)) and tissue sequence distance (D(ti)). This highly significant E(ti) - D(ti) correlation emerges when the expression divergence and protein sequence divergence are under the same tissue constraints. Second, the tissue-driven hypothesis further explains the observed significant correlation between the tissue expression distance (between the human and mouse) and the duplicate tissue distance (T(dup)) between human (or mouse) paralogous genes. In other words, between-duplicate and interspecies expression divergences covary among tissues. Third, for genes with the same expression broadness, we found that genes expressed in more stringent tissues (e.g., neurorelated) generally tend to evolve more slowly than those in more relaxed tissues (e.g., hormone-related). We conclude that tissue factors should be considered as an important component in shaping the pattern of genomic evolution and correlations.

The regulatory utilization of genetic redundancy through responsive backup circuits

Functional redundancies, generated by gene duplications, are highly widespread throughout all known genomes. One consequence of these redundancies is a tremendous increase to the robustness of organisms to mutations and other stresses. Yet, this very robustness also renders redundancy evolutionarily unstable, and it is, thus, predicted to have only a transient lifetime. In contrast, numerous reports describe instances of functional overlaps that have been conserved throughout extended evolutionary periods. More interestingly, many such backed-up genes were shown to be transcriptionally responsive to the intactness of their redundant partner and are up-regulated if the latter is mutationally inactivated. By manual inspection of the literature, we have compiled a list of such "responsive backup circuits" in a diverse list of species. Reviewing these responsive backup circuits, we extract recurring principles characterizing their regulation. We then apply modeling approaches to explore further their dynamic properties. Our results demonstrate that responsive backup circuits may function as ideal devices for filtering nongenetic noise from transcriptional pathways and obtaining regulatory precision. We thus challenge the view that such redundancies are simply leftovers of ancient duplications and suggest they are an additional component to the sophisticated machinery of cellular regulation. In this respect, we suggest that compensation for gene loss is merely a side effect of sophisticated design principles using functional redundancy.

Multiple potential germ-line helicases are components of the germ-line-specific P granules of Caenorhabditis elegans

Two components of the germ-line-specific P granules of the nematode Caenorhabditis elgans have been identified using polyclonal antibodies specific for each. Both components are putative germ-line RNA helicases (GLHs) that contain CCHC zinc fingers of the type found in the RNA-binding nucleocapsid proteins of retroviruses. The predicted GLH-1 protein has four CCHC fingers; GLH-2 has six. Both GLH proteins localize in the P granules at all stage of germ-line development. However, the two glh genes display different patterns of RNA and protein accumulation in the germ lines of hermaphrodites and males. Injection of antisense glh-1 or glh-2 RNA into wild-type worms causes some offspring to develop into sterile adults, suggesting that either or both genes are required for normal germ-line development. As these very similar glh genes physically map within several hundred kilobases of one another, it seems likely that they represent a fairly recent gene duplication event.