The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes.

BACKGROUND: Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood. METHODS: In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus. RESULTS: We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes. CONCLUSION: The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

ZmWAK02 encoding an RD-WAK protein confers maize resistance against gray leaf spot.

Gray leaf spot (GLS) caused by Cercospora zeina or C. zeae-maydis is a major maize disease throughout the world. Although more than 100 QTLs resistant against GLS have been identified, very few of them have been cloned. Here, we identified a major resistance QTL against GLS, qRglsSB, explaining 58.42% phenotypic variation in SB12×SA101 BC1 F1 population. By fine-mapping, it was narrowed down into a 928 kb region. By using transgenic lines, mutants and complementation lines, it was confirmed that the ZmWAK02 gene, encoding an RD wall-associated kinase, is the responsible gene in qRglsSB resistant against GLS. The introgression of the ZmWAK02 gene into hybrid lines significantly improves their grain yield in the presence of GLS pressure and does not reduce their grain yield in the absence of GLS. In summary, we cloned a gene, ZmWAK02, conferring large effect of GLS resistance and confirmed its great value in maize breeding.

Pushing the limits of HiFi assemblies reveals centromere diversity between two Arabidopsis thaliana genomes.

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.

Multiple independent recombinations led to hermaphroditism in grapevine.

Hermaphroditic (perfect) flowers were a key trait in grapevine domestication, enabling a drastic increase in yields due to the efficiency of self-pollination in the domesticated grapevine (Vitis vinifera L. ssp. vinifera). In contrast, all extant wild Vitis species are dioecious, each plant having only male or female flowers. In this study, we identified the male (M) and female (f) haplotypes of the sex-determining region (SDR) in the wild grapevine species V. cinerea and confirmed the boundaries of the SDR. We also demonstrated that the SDR and its boundaries are precisely conserved across the Vitis genus using shotgun resequencing data of 556 wild and domesticated accessions from North America, East Asia, and Europe. A high linkage disequilibrium was found at the SDR in all wild grape species, while different recombination signatures were observed along the hermaphrodite (H) haplotype of 363 cultivated accessions, revealing two distinct H haplotypes, named H1 and H2. To further examine the H2 haplotype, we sequenced the genome of two grapevine cultivars, 'Riesling' and 'Chardonnay'. By reconstructing the first two H2 haplotypes, we estimated the divergence time between H1 and H2 haplotypes at ∼6 million years ago, which predates the domestication of grapevine (∼8,000 y ago). Our findings emphasize the important role of recombination suppression in maintaining dioecy in wild grape species and lend additional support to the hypothesis that at least two independent recombination events led to the reversion to hermaphroditism in grapevine.

New insights into immune genes and other expanded gene families of the house fly, Musca domestica, from an improved whole genome sequence.

The house fly, Musca domestica, is a pest of livestock, transmits pathogens of human diseases, and is a model organism in multiple biological research areas. The first house fly genome assembly was published in 2014 and has been of tremendous use to the community of house fly biologists, but that genome is discontiguous and incomplete by contemporary standards. To improve the house fly reference genome, we sequenced, assembled, and annotated the house fly genome using improved techniques and technologies that were not available at the time of the original genome sequencing project. The new genome assembly is substantially more contiguous and complete than the previous genome. The new genome assembly has a scaffold N50 of 12.46 Mb, which is a 50-fold improvement over the previous assembly. In addition, the new genome assembly is within 1% of the estimated genome size based on flow cytometry, whereas the previous assembly was missing nearly one-third of the predicted genome sequence. The improved genome assembly has much more contiguous scaffolds containing large gene families. To provide an example of the benefit of the new genome, we used it to investigate tandemly arrayed immune gene families. The new contiguous assembly of these loci provides a clearer picture of the regulation of the expression of immune genes, and it leads to new insights into the selection pressures that shape their evolution.

Loss of pod strings in common bean is associated with gene duplication, retrotransposon insertion and overexpression of PvIND.

Fruit development has been central in the evolution and domestication of flowering plants. In common bean (Phaseolus vulgaris), the principal global grain legume staple, two main production categories are distinguished by fibre deposition in pods: dry beans, with fibrous, stringy pods; and stringless snap/green beans, with reduced fibre deposition, which frequently revert to the ancestral stringy state. Here, we identify genetic and developmental patterns associated with pod fibre deposition. Transcriptional, anatomical, epigenetic and genetic regulation of pod strings were explored through RNA-seq, RT-qPCR, fluorescence microscopy, bisulfite sequencing and whole-genome sequencing. Overexpression of the INDEHISCENT ('PvIND') orthologue was observed in stringless types compared with isogenic stringy lines, associated with overspecification of weak dehiscence-zone cells throughout the pod vascular sheath. No differences in DNA methylation were correlated with this phenotype. Nonstringy varieties showed a tandemly direct duplicated PvIND and a Ty1-copia retrotransposon inserted between the two repeats. These sequence features are lost during pod reversion and are predictive of pod phenotype in diverse materials, supporting their role in PvIND overexpression and reversible string phenotype. Our results give insight into reversible gain-of-function mutations and possible genetic solutions to the reversion problem, of considerable economic value for green bean production.

Mapping the genomic landscape of CRISPR-Cas9 cleavage.

RNA-guided CRISPR-Cas9 endonucleases are widely used for genome engineering, but our understanding of Cas9 specificity remains incomplete. Here, we developed a biochemical method (SITE-Seq), using Cas9 programmed with single-guide RNAs (sgRNAs), to identify the sequence of cut sites within genomic DNA. Cells edited with the same Cas9-sgRNA complexes are then assayed for mutations at each cut site using amplicon sequencing. We used SITE-Seq to examine Cas9 specificity with sgRNAs targeting the human genome. The number of sites identified depended on sgRNA sequence and nuclease concentration. Sites identified at lower concentrations showed a higher propensity for off-target mutations in cells. The list of off-target sites showing activity in cells was influenced by sgRNP delivery, cell type and duration of exposure to the nuclease. Collectively, our results underscore the utility of combining comprehensive biochemical identification of off-target sites with independent cell-based measurements of activity at those sites when assessing nuclease activity and specificity.

A maize glutaredoxin gene, abphyl2, regulates shoot meristem size and phyllotaxy.

Phyllotaxy describes the geometric arrangement of leaves and is important for plant productivity. Auxin is well known to regulate phyllotactic patterns via PIN1-dependent auxin polar transport, and studies of maize (Zea mays) aberrant phyllotaxy1 (abph1) mutants suggest the importance of auxin and cytokinin signaling for control of phyllotaxy. However, whether additional regulators control these patterns is poorly understood. Here, we report a new dominant maize mutant, Aberrant phyllotaxy2 (Abph2), in which the shoot meristems are enlarged and the phyllotactic pattern switches from alternate to decussate. Map-based cloning revealed that the Abph2 mutation was caused by transposition of a glutaredoxin gene, MALE STERILE CONVERTED ANTHER1 (MSCA1), which gained an altered expression pattern in Abph2 mutant embryos. msca1 loss-of-function mutants have reduced meristem size and revealed a novel function of glutaredoxins in meristem growth. In addition, MSCA1 interacts with a TGA transcription factor, FASCIATED EAR4, suggesting a novel regulatory module for regulating shoot meristem size.

The maize methylome influences mRNA splice sites and reveals widespread paramutation-like switches guided by small RNA.

The maize genome, with its large complement of transposons and repeats, is a paradigm for the study of epigenetic mechanisms such as paramutation and imprinting. Here, we present the genome-wide map of cytosine methylation for two maize inbred lines, B73 and Mo17. CG (65%) and CHG (50%) methylation (where H = A, C, or T) is highest in transposons, while CHH (5%) methylation is likely guided by 24-nt, but not 21-nt, small interfering RNAs (siRNAs). Correlations with methylation patterns suggest that CG methylation in exons (8%) may deter insertion of Mutator transposon insertion, while CHG methylation at splice acceptor sites may inhibit RNA splicing. Using the methylation map as a guide, we used low-coverage sequencing to show that parental methylation differences are inherited by recombinant inbred lines. However, frequent methylation switches, guided by siRNA, persist for up to eight generations, suggesting that epigenetic inheritance resembling paramutation is much more common than previously supposed. The methylation map will provide an invaluable resource for epigenetic studies in maize.

A rare SNP mutation in Brachytic2 moderately reduces plant height and increases yield potential in maize.

Plant height has long been an important agronomic trait in maize breeding. Many plant height QTLs have been reported, but few of these have been cloned. In this study, a major plant height QTL, qph1, was mapped to a 1.6kb interval in Brachytic2 (Br2) coding sequence on maize chromosome 1. A naturally occurring rare SNP in qph1, which resulted in an amino acid substitution, was validated as the causative mutation. QPH1 protein is located in the plasma membrane and polar auxin transport is impaired in the short near-isogenic line RIL88(qph1). Allelism testing showed that the SNP variant in qph1 reduces longitudinal cell number and decreases plant height by 20% in RIL88(qph1) compared to RIL88(QPH1), and is milder than known br2 mutant alleles. The effect of qph1 on plant height is significant and has no or a slight influence on yield in four F2 backgrounds and in six pairs of single-cross hybrids. Moreover, qph1 could reduce plant height when heterozygous, allowing it to be easily employed in maize breeding. Thus, a less-severe allele of a known dwarf mutant explains part of the quantitative variation for plant height and has great potential in maize improvement.

Eighteen New Candidate Effectors of the Phytonematode Heterodera glycines Produced Specifically in the Secretory Esophageal Gland Cells During Parasitism.

Heterodera glycines, the soybean cyst nematode, is the number one pathogen of soybean (Glycine max). This nematode infects soybean roots and forms an elaborate feeding site in the vascular cylinder. H. glycines produces an arsenal of effector proteins in the secretory esophageal gland cells. More than 60 H. glycines candidate effectors were identified in previous gland-cell-mining projects. However, it is likely that additional candidate effectors remained unidentified. With the goal of identifying remaining H. glycines candidate effectors, we constructed and sequenced a large gland cell cDNA library resulting in 11,814 expressed sequence tags. After bioinformatic filtering for candidate effectors using a number of criteria, in situ hybridizations were performed in H. glycines whole-mount specimens to identify candidate effectors whose mRNA exclusively accumulated in the esophageal gland cells, which is a hallmark of many nematode effectors. This approach resulted in the identification of 18 new H. glycines esophageal gland-cell-specific candidate effectors. Of these candidate effectors, 11 sequences were pioneers without similarities to known proteins while 7 sequences had similarities to functionally annotated proteins in databases. These putative homologies provided the bases for the development of hypotheses about potential functions in the parasitism process.

tassel-less1 encodes a boron channel protein required for inflorescence development in maize.

tassel-less1 (tls1) is a classical maize (Zea mays) inflorescence mutant. Homozygous mutant plants have no tassels or very small tassels, and ear development is also impaired. Using a positional cloning approach, ZmNIP3;1 (a NOD26-like intrinsic protein) was identified as the candidate gene for tls1. The ZmNIP3;1 gene is completely deleted in the tls1 mutant genome. Two Mutator-insertional TUSC alleles of ZmNIP3;1 exhibited tls1-like phenotypes, and allelism tests confirmed that the tls1 gene encodes ZmNIP3;1. Transgenic plants with an RNA interference (RNAi) construct to down-regulate ZmNIP3;1 also showed tls1-like phenotypes, further demonstrating that TLS1 is ZmNIP3;1. Sequence analysis suggests that ZmNIP3;1 is a boron channel protein. Foliar application of boron could rescue the tls1 phenotypes and restore the normal tassel and ear development. Gene expression analysis indicated that in comparison with that of the wild type or tls1 plants treated with boron, the transition from the vegetative to reproductive phase or the development of the floral meristem is impaired in the shoot apical meristem of the tls1 mutant plants. It is concluded that the tls1 mutant phenotypes are caused by impaired boron transport, and boron is essential for inflorescence development in maize.

Current status of community resources and priorities for weed genomics research.

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility.

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

The gyrfalcon (Falco rusticolus) genome.

High-quality genome assemblies are characterized by high-sequence contiguity, completeness, and a low error rate, thus providing the basis for a wide array of studies focusing on natural species ecology, conservation, evolution, and population genomics. To provide this valuable resource for conservation projects and comparative genomics studies on gyrfalcon (Falco rusticolus), we sequenced and assembled the genome of this species using third-generation sequencing strategies and optical maps. Here, we describe a highly contiguous and complete genome assembly comprising 20 scaffolds and 13 contigs with a total size of 1.193 Gbp, including 8,064 complete Benchmarking Universal Single-Copy Orthologs (BUSCOs) of the total 8,338 BUSCO groups present in the library aves_odb10. Of these BUSCO genes, 96.7% were complete, 96.1% were present as a single copy, and 0.6% were duplicated. Furthermore, 0.8% of BUSCO genes were fragmented and 2.5% (210) were missing. A de novo search for transposable elements (TEs) identified 5,716 TEs that masked 7.61% of the F. rusticolus genome assembly when combined with publicly available TE collections. Long interspersed nuclear elements, in particular, the element Chicken-repeat 1 (CR1), were the most abundant TEs in the F. rusticolus genome. A de novo first-pass gene annotation was performed using 293,349 PacBio Iso-Seq transcripts and 496,195 transcripts derived from the assembly of 42,429,525 Illumina PE RNA-seq reads. In all, 19,602 putative genes, of which 59.31% were functionally characterized and associated with Gene Ontology terms, were annotated. A comparison of the gyrfalcon genome assembly with the publicly available assemblies of the domestic chicken (Gallus gallus), zebra finch (Taeniopygia guttata), and hummingbird (Calypte anna) revealed several genome rearrangements. In particular, nine putative chromosome fusions were identified in the gyrfalcon genome assembly compared with those in the G. gallus genome assembly. This genome assembly, its annotation for TEs and genes, and the comparative analyses presented, complement and strength the base of high-quality genome assemblies and associated resources available for comparative studies focusing on the evolution, ecology, and conservation of Aves.

The northern corn leaf blight resistance gene Ht1 encodes an nucleotide-binding, leucine-rich repeat immune receptor.

Northern corn leaf blight, caused by the fungal pathogen Exserohilum turcicum, is a major disease of maize. The first major locus conferring resistance to E. turcicum race 0, Ht1, was identified over 50 years ago, but the underlying gene has remained unknown. We employed a map-based cloning strategy to identify the Ht1 causal gene, which was found to be a coiled-coil nucleotide-binding, leucine-rich repeat (NLR) gene, which we named PH4GP-Ht1. Transgenic testing confirmed that introducing the native PH4GP-Ht1 sequence to a susceptible maize variety resulted in resistance to E. turcicum race 0. A survey of the maize nested association mapping genomes revealed that susceptible Ht1 alleles had very low to no expression of the gene. Overexpression of the susceptible B73 allele, however, did not result in resistant plants, indicating that sequence variations may underlie the difference between resistant and susceptible phenotypes. Modelling of the PH4GP-Ht1 protein indicated that it has structural homology to the Arabidopsis NLR resistance gene ZAR1, and probably forms a similar homopentamer structure following activation. RNA sequencing data from an infection time course revealed that 1 week after inoculation there was a threefold reduction in fungal biomass in the PH4GP-Ht1 transgenic plants compared to wild-type plants. Furthermore, PH4GP-Ht1 transgenics had significantly more inoculation-responsive differentially expressed genes than wild-type plants, with enrichment seen in genes associated with both defence and photosynthesis. These results demonstrate that the NLR PH4GP-Ht1 is the causal gene underlying Ht1, which represents a different mode of action compared to the previously reported wall-associated kinase northern corn leaf blight resistance gene Htn1/Ht2/Ht3.

Near-gapless genome assemblies of Williams 82 and Lee cultivars for accelerating global soybean research.

Complete, gapless telomere-to-telomere chromosome assemblies are a prerequisite for comprehensively investigating the architecture of complex regions, like centromeres or telomeres and removing uncertainties in the order, spacing, and orientation of genes. Using complementary genomics technologies and assembly algorithms, we developed highly contiguous, nearly gapless, genome assemblies for two economically important soybean [Glycine max (L.) Merr] cultivars (Williams 82 and Lee). The centromeres were distinctly annotated on all the chromosomes of both assemblies. We further found that the canonical telomeric repeats were present at the telomeres of all chromosomes of both Williams 82 and Lee genomes. A total of 10 chromosomes in Williams 82 and eight in Lee were entirely reconstructed in single contigs without any gap. Using the combination of ab initio prediction, protein homology, and transcriptome evidence, we identified 58,287 and 56,725 protein-coding genes in Williams 82 and Lee, respectively. The genome assemblies and annotations will serve as a valuable resource for studying soybean genomics and genetics and accelerating soybean improvement.

A happy accident: a novel turfgrass reference genome.

Poa pratensis, commonly known as Kentucky bluegrass, is a popular cool-season grass species used as turf in lawns and recreation areas globally. Despite its substantial economic value, a reference genome had not previously been assembled due to the genome's relatively large size and biological complexity that includes apomixis, polyploidy, and interspecific hybridization. We report here a fortuitous de novo assembly and annotation of a P. pratensis genome. Instead of sequencing the genome of a C4 grass, we accidentally sampled and sequenced tissue from a weedy P. pratensis whose stolon was intertwined with that of the C4 grass. The draft assembly consists of 6.09 Gbp with an N50 scaffold length of 65.1 Mbp, and a total of 118 scaffolds, generated using PacBio long reads and Bionano optical map technology. We annotated 256K gene models and found 58% of the genome to be composed of transposable elements. To demonstrate the applicability of the reference genome, we evaluated population structure and estimated genetic diversity in P. pratensis collected from three North American prairies, two in Manitoba, Canada and one in Colorado, USA. Our results support previous studies that found high genetic diversity and population structure within the species. The reference genome and annotation will be an important resource for turfgrass breeding and study of bluegrasses.

Improved pearl millet genomes representing the global heterotic pool offer a framework for molecular breeding applications.

High-quality reference genome assemblies, representative of global heterotic patterns, offer an ideal platform to accurately characterize and utilize genetic variation in the primary gene pool of hybrid crops. Here we report three platinum grade de-novo, near gap-free, chromosome-level reference genome assemblies from the active breeding germplasm in pearl millet with a high degree of contiguity, completeness, and accuracy. An improved Tift genome (Tift23D2B1-P1-P5) assembly has a contig N50 ~ 7,000-fold (126 Mb) compared to the previous version and better alignment in centromeric regions. Comparative genome analyses of these three lines clearly demonstrate a high level of collinearity and multiple structural variations, including inversions greater than 1 Mb. Differential genes in improved Tift genome are enriched for serine O-acetyltransferase and glycerol-3-phosphate metabolic process which play an important role in improving the nutritional quality of seed protein and disease resistance in plants, respectively. Multiple marker-trait associations are identified for a range of agronomic traits, including grain yield through genome-wide association study. Improved genome assemblies and marker resources developed in this study provide a comprehensive framework/platform for future applications such as marker-assisted selection of mono/oligogenic traits as well as whole-genome prediction and haplotype-based breeding of complex traits.